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Abstract. For any i? > 0, we construct a sequence of bounded 
potentials V^f , n G N, supported in an annular region Bout\Binn C 
M'^, which act as approximate cloaks for solutions of Schrodinger's 
equation at energy E: For any potential Vq G L°°{Binn) such that 
E is not a Neumann eigenvalue of — A + Vq in -Bm„, the scatter- 
ing amplitudes ay^+yE (i?, 0, w) as n ^ oo. The thus 
not only form a family of approximately transparent potentials, 
but also function as approximate invisibility cloaks in quantum 
mechanics. On the other hand, for E close to interior eigenval- 
ues, resonances develop and there exist almost trapped states con- 
centrated in Binn- We derive the from singular, anisotropic 
transformation optics-based cloaks by a de-anisotropization pro- 
cedure, which we call isotropic transformation optics. This tech- 
nique uses truncation, inverse homogenization and spectral theory 
to produce nonsingular, isotropic approximate cloaks. As an in- 
termediate step, we also obtain approximate cloaking for a general 
class of equations including the acoustic equation. 



1. Introduction 

A fundamental problem is to describe the scattering of waves by a 
potential, as governed by the time-independent Schrodinger equation 
at energy E > 0, 

(1) {-A + V{x))^{x) = Eij{x), xgM^ 

ipix) = exp{iE^^'^x- 6) + ipsc{x), 

where 6 G §^~^ and ipsc{x) satisfies the Sommerfeld radiation condition. 
We restrict ourselves in this paper to compactly supported potentials 
V, so 



ipscix) = CdE 4 exp{iE''^\x\) + o[\x\ 2 

I T I 5 \ / 



as \x\ 00. 
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The function av{E,u!,6) is the scattering amplitude at energy E of 
the potential V. The associated inverse scattering problem consists of 
trying to determine V from the scattering amplitude, or measurements 
of waves at the boundary of some region Q containing the support of 
V. 

Recently, Zhang, et ai, ^65j described quantum mechanical cloaking 
at any fixed energy E. Their construction starts with a homogeneous, 
isotropic mass tensor and a potential V = 0, and subjects this pair to a 
singular change of variables. This "blowing up a point" transformation 
had been used in EH] to produce conductivities that hide objects 
from detection by electrostatic measurement^, and was subsequently 
used to describe the same phenomenon for electromagnetic waves |54j : 
one now refers to a specification of material parameters having this 
effect as a cloak. The cloaking Schrodinger equation in [65], which has 
an anisotropic, singular mass tensor, is equivalent with the Helmholtz 
equation (at frequency u = \fE) for an associated singular Riemannian 
metric, and thus covered by cloaking for the Helmholtz equation in 3D, 
as we analyzed in [20l Sec. 3]. Similarly, cloaks for acoustics in 3D have 
been described in (TUl fTB|^ : again, these are in fact direct consequences 
of cloaking for the Helmholtz equation in 3D, cf. [24J. 

Thus, for ideal 3D cloaking in each of scalar optics, quantum mechanics 
and acoustics, one knows from PU] that any finite energy distributional 
solution decouples into a sum of a wave on the exterior of the cloak, 
unaffected (in terms of scattering or boundary measurements) by the 
cloak, and a wave within the cloaked region satisfying the Neumann 
boundary condition at the cloaking surface. Hence, if E is not a Neu- 
mann eigenvalue, then the wave must vanish within the cloaked region 
and cloaking works as advertised. On the other hand, if E is an eigen- 
value, the cloaked region supports interior resonances, or trapped states. 
This is an unphysical situation, since the Dirichlet problem on f2 no 
longer has unique solutions, and this can be considered as a failure of 
cloaking, per se. However, what emerges from this failure in the setting 
of approximate cloaking described below is, we believe, quite useful. 

In this paper, we construct, for each energy E, a family {V^^jJ^^^ of po- 
tentials, supported in an annulus Bout \ -Bmn C M^, which are not only 
almost transparent, in that the scattering amplitudes ayE{E, -,•)—> 
as n ^ oo, but also act as approximate cloaks for potentials supported 



The 2D version has the same electrostatic cloaking property j36| . 
^See [13] for the 2D case, and [3S1 for more regarding elastic and acoustic 
cloaking. 
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in the inner ball Binn'- for any potential Vq G L°° (Binn) for which E 
is not an eigenvalue of —A + Vq in Binn, the scattering amplitudes 
avQ+v^{E, ■, ■) — > as n ^ oo, as well. There are also approximate 
versions of the interior resonances supported by the ideal cloak: there 
exist energies close to the Neumann eigenvalues of the cloaked region 
Binn, near which there are waves largely concentrated in Binn, which 
we call almost trapped states. Furthermore, the quality of the approxi- 
mate cloaking degrades as we move towards these energies, with waves 
being able to penetrate -Bm„, and the presence of this region and Vq 
detectable by scattering or boundary measurements. In addition, we 
are able to include a magnetic potential in the Schrodinger equation; 
this allows one to switch between the approximate cloak and almost 
trapped state alternatives by application of a suitable homogeneous 
magnetic field. We have given physical applications of this elsewhere 
[251 [26]. 

To put these results in a mathematical context, recall that in several 
dimensions numerous uniqueness results are available for the inverse 
scattering and boundary value problems, subject to some regularity 
assumptions, in terms of the scattering data at a fixed energy E and 
the closely related Dirichlet-to-Neumann (DN) operator. Physically, 
these correspond to far- and near-field measurements, resp. In fact, 
for compactly supported potentials the scattering data is equivalent to 
the DN operator measured on the boundary of a domain containing 
the support of the potential [6]. In dimensions d > 3, Sylvester and 
Uhlmann [62] showed that a smooth potential is determined by the 
associated DN map, and this was extended to L°° in [49j. Lavine 
and Nachman [52] and Chanillo [S] generalized this further, showing 
that, if a potential V G L2(M'^) or has small norm in the Fefferman- 
Phong space Fp, p > resp., then the DN map determines V. For 
d = 2, Nachman [48j proved that uniqueness holds if V has the special 
form V = —A{y^)/y/^ arising from the gauge transformation of the 
conductivity equation to the Schrodinger equation, with 7 G W'^'^, p > 
1, while Sun and Uhlmann [61] showed that uniqueness holds for a 
generic class of potentials. Very recently, Bukhgeim[7] has shown that 
uniqueness holds for all potentials in L^Q^p(M^), p > 2. 

On the other hand, for each positive energy E, Regge [57j constructed 
(noncompactly supported) potentials which are transparent, i.e., for 
which the scattering amplitude av{E,9,uj) = ao{E,9,uj) = 0. See also 
[5ni[59l|30], where the last gives a construction of transparent potentials 
in the Schwartz class iS(M^). In the physics literature, Hendi, Henn and 
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Leonhardt [31] have described central, i.e., radial, potentials which are 
transparent on the level of the ray geometry. 

We construct the families {V^} of approximately cloaking potentials 
by means of a result of independent interest. The use of changes of vari- 
ables to produce novel optical effects on waves or to facilitate computa- 
tions has been considered in the physics literature, (se e|, e.g., Dolin [17] 
or more recently Ward and Pendry [63j), and is now generally referred 
to as transformation optics (TO). However, to produce cloaking and 
other extreme effects, nonsingular changes of variables are insufficient. 
The recently proposed plans for cloaking are based on singular transfor- 
mations and consist of medium parameters which are both anisotropic 
and singula^ whether for the conductivity (electrostatics) [281 ESj; 
index of refraction (Helmholtz) [HI [201 ISS] ; permittivity and perme- 
ability (Maxwell) [SI [20]; mass density (acoustic) [H [TOl [151 [2l] ; or 
effective mass (Schrodinger) [65] . Physical realization of such designs is 
now potentially feasible due to the rapidly developing area of metama- 
terials, but the singularity and extreme anisotropy make characterizing 
and fabricating the materials to implement such designs an enormous 
obstacle to manufacturing invisibility devices. These same remarks are 
valid for other TO designs, such as [1^1 [HI HH [211 [2SIISS| , some of which 
are singular. 

We propose here a general method, which we refer to by the oxymoronic 
isotropic transformation optics, for dealing with both the anisotropy 
and singularity of TO material parameters. We describe this in detail 
in the context of cloaking, but it should be applicable to a wider range 
of TO designs. We in fact derive the quantum mechanical approximate 
cloaks from approximate cloaks for a general class of equations that 
includes the acoustic equation. Using ideas from Nguetseng [JT] , Allaire 
[T] , Cherkaev [13] and elsewhere, we show how to find cloaking material 
parameters that are at once both isotropic and nonsingular, at the 
price of replacing perfect (ideal) cloaking with approximate cloaking 
(of arbitrary accuracy). This method, starting with transformation 
optics-based designs and constructing approximations to them, first 
by nonsingular, but still anisotropic, material parameters, and then 
by nonsingular isotropic parameters, appears to be a very flexible tool 
for creating physically realistic theoretical designs, easier to implement 
than the ideal ones due to the relatively tame nature of the materials 



We thank A. Kildishev for this reference. 

^By singular, we mean that at least one of the eigenvalues goes to zero or infinity 
at some points. 
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needed, yet (up to an arbitrarily small error) essentially capturing the 
desired effects on wave propagation. 

We review the ideal electrostatic cloak of [2H], and extend this in Sec. 
12.11 to a class of equations which will allow us to deal with both gen- 
eral acoustic and magnetic Schrodinger equations. These equations 
have coefficients which, borrowing the terminology from acoustics, we 
refer to as mass density and bulk modulus and which are singular at 
the cloaking surface S, i.e., the interface between the cloaked and un- 
cloaked regions. The desingularization process begins in Sec. 12.21 by 
truncating the mass density away from E, on the outer side of which 
the ideal cloaking parameters are singular. (Similar truncations have 
been considered before in the context of cylindrical or 2D cloaking, cf. 
[581 [221 [TH [36].) We show in Sec. [23] that the Dirichlet forms for the 
ideal cloaks can be well-approximated (in the sense of F- convergence) 
by these truncations. We then desingularize the bulk modulus in Sec. 
[M 

So far, the approximately cloaking mass densities are still anisotropic. 
As is well known in effective medium theory, homogenization of iso- 
tropic material parameters may lead to anisotropic ones [45]; in Sees. 
[3]and[ll we use this phenomenon in reverse, showing that the Dirichlet 
forms obtained in Sec. [2] can be well-approximated by those for certain 
nonsingular isotropic conductivities, which thus provide approximate 
cloaks for the general class of acoustic-like equations. These then allow 
us to obtain in Thm. 15.11 families of approximate quantum cloaks. In 
Sec. [6] we study failure of cloaking near exceptional energies, mirroring 
the failure of the existence of unique solutions for the ideal cloak at 
Neumann eigenvalues of the cloaked region; see the Remark at the end 
of the section. Further physical applications, including a new type of 
ion trap, can be found in Finally, numerical simulations are 

presented in Sec. [7] 

AG and GU were supported by the US NSF, YK by EPSRC, ML by 
Academy of Finland (CoE project 213476), and GU by a Walker Family 
Professorship. The authors thank A. Cherkaev and V. Smyshlyaev for 
useful discussions on homogenization. 

2. Approximating cloaking material parameters by 
nonsingular anisotropic parameters 

We will denote by B{a,R) the ball of radius R centered at a in M^, 
sometimes denoted simply B{R) when centered at the origin, O. Let 
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Ml = 5(0i,3) C and Ms = B{02,1), with O^Os being two 
copies of O, considered as disjoint compact manifolds with boundary; 
set M = Ml U Ms. Also, let = B{3) C M^ and : Mi \ {0;} ^ 
n\^{l) C be the map 



Together, these form a surjective map F = [F^^F"^) from the cloaking 
manifold (or virtual space) M \ { } to the cloaking device (or physical 
space) Vt. By a conductivity we mean a measurable map with values in 
the symmetric non-negative M?^'^ matrices. Let 70 = 1 be the constant 
isotropic conductivity on and define the conductivity cti on VL as 



(4) ai = F,So, xen\B{l), ai = 2^270, xeB{l), 



which has a singularity on the cloaking surface S := (9-8(1), both in that 
one of the eigenvalues (corresponding to the radial direction) tends to 
as r \ 1 and that there is a jump discontinuity across the sphere S. 
This conductivity cxi is, up to the radius of f2 and the factor 2 in the sec- 
ond of formulae (jlj), used here for technical reasons, the one introduced 
in [2HIE9] and shown to be indistinguishable from 70, vis-a-vis electro- 
static boundary measurements at dQ. In fact, cri|B(i) can be replaced 
by any smooth, non-degenerate anisotropic conductivity tensor and its 
values will be undetectable at dfl. The same construction of (Ti\n\B{i), 
applied instead to the electric permittivity e and magnetic permeability 
/i in Maxwell's equations, was proposed in |54j (see also ^1]) to cloak 
the region -B(l) from observation by electromagnetic waves at positive 
frequency; we thus refer to cxi as a cloaking conductivity, and following 
the physics literature, we will refer to {Q, o"i) as the ideal cloak. 

This gives rise to the Dirichlet problem for the singular conductivity 
equation, 

(5) V • ctiVm = 0, in fi, u\dn = h, 

and to the corresponding singular conductivity operator A, 



(2) 



F^{x)=x, for 2 < |x| < 3, 




for < Ixl < 2. 




(6) 



Au := -g'^^^V ■ (TiVm, 
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that we consider with Dirichlet boundary condition u\dQ. = 0. Here, we 
use the singular Riemannian metric {gjk)\k=i associated to the conduc- 
tivity (Ti, namely, 

(7) 9"'9^' = a{\ ^=|det[^^^]|-i = |detK]r, 

c/. |28l |29]. We denote by g both the metric and the corresponding 
scalar function, the meaning being clear from the context. 

A rigorous definition of the meaning of ([5]) and the operator A is given 
in the following sections. In particular, the operator A is self-adjoint 
on Li^giyt)^ the weighted L^-space defined using the weight g^f^. For 
a general weight w{x) > 0, we denote by L?{yL,wdx) = L?{wdx) the 
weighted space, 

L'^{VL,wdx) = {m : f2 — > C measurable, ||M||i2(^^2.-) := / \u\^wdx < oo}. 

For simplicity, we denote = L'^g{^) '■= L'^{^1, g^dx), the natural 
space for the metric g, and the norm in this space by || ■ \\g. Note that 

(8) L'^{n)cLl{n), \\u\\g < V8\\u\\. 
We also use the Sobolev spaces 

Hlin) = {ue Llin) : M|n\s G \ E), / atd.ud^dx < oo}, 

K,(^) = e Hlin) ■.v\9n = 0}. 

Here and below, we use Einstein summation convention, summing over 
indices j and k appearing both as sub- and super-indices. Observe that 
H\Q) C Hl{Q) and 

(9) = (||m||^2 + / aC'djudi^dx)^ < V8\\u\\hi- 
^ ^ Jn 

Throughout, we also use the following standard terminology: when 
considering convergence of sequences {xn}'^=i in a Hilbert space H, we 
say that x„ converges strongly to x in H, if — ^ 0, as n — oo, 
while Xn converges weekly to x, if (x„ — x, ^ 0, as n ^ oo, for any 
yeH. 

2.1. Ideal cloaking for more general equations. In this paper we 
treat equations more general than ([5]), which are important for physical 
applications, cf. [25l [26]. To this end, consider a family of equations 
that simultaneously includes both the magnetic Schrodinger equation 
and the acoustic equation; for simplicity, we will refer to equations of 
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this type as acoustic. Let q and b = (&i,&2,^3) be scalar- and vector- 
valued functions, which eventually will represent the electric and mag- 
netic potentials; we will assume that 

b e L°°(l]; M^), q G M), g > 0. 

We note that the condition g > is merely a convenience, since for gen- 
eral q G L°°(f2;M) we can always add a constant to achieve positivity, 
and of course this just shifts the spectrum. 

To deal rigorously with the elliptic boundary value problem 

(10) - ^"^/^(V + ib) • ai(V + ib)u + qu~ \u = f, u\an = K 

which has singular coefficients, we consider the corresponding quadratic 
form. In the following, we use the notation 

V^ = V + z6(x). 
Lemma 2.1. The quadratic form 

(11) ai[u,u] = / (TiV^u -Vhidx + I qg^/'^\u\'^ dx, 

Jq\t. Jn\T. 

defined in the domain T'(ai) = HQg{Q) is closed. Moreover, the em- 
bedding V^ai) > L'^{fl) is compact, and there is Cb > I so that 

(12) C;;^\\u\\li (f^) < \\u\\l + ai[u,u] < Cb\\u\\l^ 



Proof. To prove the assertion, we need to show two facts. First, we 
need to prove that ai[ti,u] < cxo for u G Hq^^Q). Second, we have to 
prove the closedness of the form ai on HQg{Q). 

We start with the fact |20] that the map F^, is unitary from {Mi) © 
L\M2) to and from H^{Mi)®H\M2) to Hlg{n), respectively. 

For V = (fi, ^2) and u = F^v, we have 

(13) {vi, V2) e L\Mi) ® L\M2) mueLlin)- 

{vi, V2) G Hl{Mi) © H\M2) iff G Hlg{9). 

Hence, as H^{Mj), j = 1,2 is compactly embedded into L'^{Mj), we 
see that the space Hg{fl) is compactly embedded into L'^{^). 
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Using definition (jlj) and the transformation rule for 1— forms, we see 
that 

(14) ai[u,u\ = a[v,v] = a^[vi,Vi] + a'^[v2,V2]; 

a^[vi,Vi]= / I Vfi + f ipdx + / Ki\vi\'^dx, 



J Ah J Ml 

.2 



a [v2,V2] = 2 / \Vv2 + iP2V2\ dx + 8 K2\v2\ dx. 

Here the 1— forms [3i = P\mi, P2 = P\m2, and functions ki = k\mi, 
1^2 = f^\M2i a^^6 given by 6 = F^[3, q = F^k, that is, 

It follows from ([2]), ([3]) that [32, 1^1 and H2 are bounded, but /5i has a 
singularity at a; = 0, of the order l/\x\, and 

(15) < 

\x\ 

Let us consider a^[f2, ^2] as an unbounded non- negative quadratic form 
in L^(M2), with domain T>{a'^) = H^{M2). Then is closed. The qua- 
dratic form a^[wi,t>i] requires further analysis. We consider a^fvi,!"!] 
as an unbounded non- negative quadratic form in L^(Mi) having the 
domain 

(16) V{a') = {v, e L\Mi) : V^;i + iPiVi e L^Mi), viIqm, = 0}. 

Note that the condition Vui -|- i/3ifi G L^(Mi) implies that the trace 
viIqm, is well-defined. By |^[60], C^(Mi) is dense on V{a^) and 
on V{a^) is a closed, non-negative quadratic form. 

By Hardy's inequality it follows from f[T^ that 

(17) IIA^^iIIl2(Mi) < C'II&IIl- II^^iIIhi(Mi), 

so that 

(18) ^o(Mi) C V{a^) is dense. 

Let vi e V{a^). Then, by [40j \vi\ G i^o(^i)- Using ([17]) again, we 
obtain 

||/?1^^i||l2(Mi) = II Al^^ll ||l2(Mi) < 00. 

Thus, by definition (fT6|) . this yields that Vfi G L^(Mi), and hence 
t;i G Hl{Mi). This shows that ©(a^) C Hl{Mi) which, together with 
(fT8|) . implies that Vi^a^) = Hq{Mi). Thus, the domain of the closed 
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form a is T>[ol) = Hq{Mi) © H^{M2). Using the transformation rule 
ffn|) for u = F^v, and f[T^ . we conclude that ai[u,u] < oo for all u G 
HQg{Q) and that the quadratic form Oi with domain V{ai) = i?dg(n) 
is closed. The inequality (fT2l) follows from this by the open mapping 
theorem. □ 

By the theory of quadratic forms [33], the closed quadratic form ai[u, u] 
defines an associated non- negative self-adjoint operator Ai in L'^i^), 
having domain 

{19)V{Ai) := {u e Hlg{n) : there is / G L2(fi) such that 

ai[u,v] = {f,v)g for v G Hlg{n)}, 
that is, for u G V{Ai) we have, for all v G Hq^^Q), 

(20) / ai\/\i-Wh^dx+ / qg^/\vdx= / g^^'^fvdx 
Jn\T. ic\E Jq\t, 

with Aiu = / G L^(r2). By Lemma [2. ![ the spectrum of Ai consists of 
discrete eigenvalues with finite dimensional eigenspaces. 

Observe that, formally integrating by parts in (1201) . we come to ( fTOl) 
with h = 0. To better understand the nature of the operator Ai, we 
give also an alternative definition in the case when b G C(i7;M'^). In 
this case, we use the fact that the map : (p ^ a{ {dk(p + ibk4>), 
defined initially for G C^{Q), has a bounded extension 

where A^(f2;M) denotes the space of Borel measures on fl, cf. [201 
Lem. 3.2]. For b G C(i7;]R^), an equivalent definition of the operator 
Ai is then 

(21) Ai{u) := -g-^/\dj + ibj)Diu + qu, 

V{Ai) = {ue Hly{n) : {dj + ibj)Dlu G L'^iVl^g-^'^dx)}. 

2.2. Approximate cloaking by truncation. For any 1 < i? < 2, 

consider the nonsingular truncations of ffTOl) . 

(22) - V''- aRV\ + qg^/'^u - \g^''^u = mVL, 

u\dn = h, 

where are measurable anisotropic conductivities in Q satisfying 

(23) cT/j I n\B(o,5/2) = 7o, ctrIbcci) = 270, 

lim aji{x) = cr(x), Ci{R - 1)70 < aR{x) < C270, for x e n, 

c^R^ix) < aR^{x), for i?2 > ^1, 
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for some Ci < 1/2, C2 > 2. For instance, we can choose 
(24) aR{x) 



F,-fQ, xen\B{R), 
270, X G B{R). 



We note that, by smoothing the conductivities fl2^ . it is possible to 
construct conductivities G (7^(0; M^) , which we use in Sec. [5l 
satisfying (1251) . 

We denote the solution of fl22]) by m = u^. Note that, for 6 = 0, g = 0, 
equation fl2^ is an acoustic equation with mass density and ^fi/A; 
modulus g^/'^ with = V^; in the quantum mechanical setting, 
corresponds to the effective mass and (g — X)g^^'^ to the potential. By 
abuse of notation, even for 6 7^ and g 7^ 0, we will refer to as the 
mass density. 

Observe that, for each i? > 1, the mass density ctr is nonsingular, i.e., 
is bounded from above and below with, however, the lower bound going 
to as \ 1. Moreover, for any a; G M^, the symmetric matrix valued 
function R ^ <Jr{x) is increasing @ as a function of i?, and therefore 
decreases as -R \ 1. Nonsingular regularizations or truncations of 
singular ideal cloaks have previously been considered in [581 [221 [36] . 

To motivate the treatment here, consider for R > 1 the Dirichlet-to- 
Neumann (DN) map : H^/'^{dn) H'^^'^^dQ) that maps 

(25) : u\dn ^ (d^u + iu ■ bu)\dn, 

where u solves f[22l) . The DN map corresponds to the Dirichlet-to- 
Neumann quadratic form, which by abuse of notation is denoted also 
by Ai 

(26) A'^[h] = [ {A'^h){x)h{^dS{x), 

J an 

where we denote h = u\qq G H^^'^{dVL)] for A < 0, the Dirichlet-to- 
Neumann form may be also represented as 

(27) A\\h]=M{aR[u,u]-X\\u\Q, 

where infimum is taken over all u G H^{VL) with u\qq, = h. However, 
to treat general A, and the general class of equations (!22l) . we will use 
the definition f[3Tl) below. 



^Note that, due to the behavior of the eigenvalues of ai , this simple but important 
monotonicity property fails in the 2D case, which we are not treating. 
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Returning to equation ( l22l) . note that for A < and R > 1 the solu- 
tion can be obtained from the minimization problem for the quadratic 
functional associated to the sesquilinear functional 

(28) aif [u,v]= / aii{x)'V''u{x)-'V^v{x) dx + / q{x)g^^'^{x)u{x)v{x) dx 
Jq Jq 



Moreover, we have 



u 



R 



argmin (aR[M,M] - A||u||^) 



where minimization is taken over u G H^{Q) such that u\dQ = h. 
Observe that the DN form (1271) is related to qr defined by (1281) . 

Thus, the solution of the acoustic equation ( l22l) with the Dirichlet 
boundary value h is the minimizer of aR[u,u] — A||m||^. 

Now consider the quadratic form aR[u,u] restricted to n G HqIQ). 

Lemma 2.2. For R > 1, the sesquilinear form with domain 
T^{.o,r) = Hq{Q), is closed and positive definite on L'^{^)- 

Proof. Clearly an is non-negative. The fact that is closed on Vr = 
Hq{Q) follows from the same considerations as those in Lemma l2.ll 
Next we show that V^u = iff m = which, due to the compactness 
of the embedding Hq{Q) ^ L^, implies the positive-definiteness of a^. 
Assume, on the contrary, that there is a non-zero ip G Hq{Q) such that 
V'^tp = in Q. Continue b and by to M'^ \ Q. and also extend 
to \ as 70 = 1 in \ b{3). Then ip G H^{R^) and VV = in 
M^ and hence V'' ■ VV = in \ {O}. Using unique continuation 
for elliptic equations with non-smooth first order terms, see [ll[3l], we 
deduce that ?/' = OinM3\{0}. □ 

As (Tr is bounded from below. Lemma 12.21 implies that there is cr> 
such that 

(29) / cxrV% ■ Vhidx > CR||M||^2(f^), u G H^i^). 
Jn 

Similarly to ( |T9|) . (120|) . the sesquilinear form qr defines a self-adjoint 
operator Ar in L^, and, analogously to ( 12T|) . we have 

(30) Ar{u) := -g-'/^V' ■ (aRV'u) + qu, 

V{Ar) = {ue H^{n) : V' ■ (aRV'u) G L^{n, g-^/^dx)}. 
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Using the operator Ar, i? > 1 we see that, for A ^ spec(.4_R), 

(31) A^nh = {d, + iv-h) (mo + ^^r) 

:= 7^R(A) (V'' ■ (ctrV'uo) - quo + \uo) . 

Here uq G H^{^1) with supp(Mo) C ^l\B{0, 2) satisfies UqIqq = h, while 
71r{X) is the resolvent, 

(32) -JZniX) = {Ar - XI)-' : H-\Q) ^ L^Q). 

We will use the notation fl52]) for i? = 1 as well as i? > 1, in which case 
7^l(A) : H-'{n) Ll{n). Note that the right-hand side in (EI]) should 
be understood in the sense of the pairing if~^/^((9fi) x H^^^'^(d^l), 

/ {{d^ + iu ■ b){uo + v^))ip dS = aR[uo + v^, tp] ~ X{uo + v^, tp) , 
Jan 

for all ip G H'{Q), and we again consider qr on the whole of H'{Q). 

Next, consider the DN map in the ideal case R = 1. Following [20] . 
we say that m is a finite energy solution of the boundary value problem 
([ID]) with h G H^^^dQ) and / G H-\n), Hue i/](fi) and 

(33) a,[u-uo,(l)] = - I {<yiV'uo-V^+{qg''^uo-Xg"^u- f)^dx, 

Jn 

for every (f) G H^g^Q). Here Hg^{Q) is the dual space to Hq^^Q). 
On the other hand, 

(34) -{V + ipi)-{V + ipi)vi + K,vi-Xv^ = f,, in Ml, 

VlldMi = h, 

(35) -^(V + 1(32) ■ (y + iP2)v2 + K2V2 - Xv2 = h, inMs, 

(9,, + il^- P2)v2\dM2 = 0, 

are satisfied in the weak sense if 

(36) a[v — -uo, -^A] = — / (cj V^t;o ■ V^'?/' + (kjUo — Xv — fj)ip) dx, 

for all ifj G i^o(^i) © H\M2), where Ci = 1, C2 = |. Here G 
if^(Mi) © H^{M2) is supported in Mi \ 5(2) and satisfies Vo\dMi = h. 

Lemma 2.3. Let h G i7^/^((9fi) and f G Hg'^{n). A function u G 
Hg{Q) is a finite energy solution in the sense (E^j of the boundary 
value problem / flOj) if and only if v = (fi,f2) G H^{Mi) © H^{M2), 
u = F^:V satisfies equations (34^ and (E^ m i/ie tfea/i; sense (E^j wt/i 

/ = i^*(7i,72). 
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In particular, for any A G M the Cauchy data, on dfl, of solutions to 
^M) satisfy 

(37) {{u\an, [d^ + iv- h)u\dQ) :-V*- aiV^u + qg^'^u = g^'^Xu m Vl} 
= {{vi\dn,{du + ii^- Pi)vi\dn) ■■-V'^^-V'^'^vi + KiVi = Xvi in n} 

Proof. By Lemma 12.11 and formulae (fT3l) and (IT^ we see that u G 
V{ai) if and only if u = F^v, v = {vi, V2) G V^a), and 

(38) ai[^,0] = for 0, G P(a), = ^ = F.V'- 



Let vq G H\Mi) © H^{M2) be supported in Mi \ fi(2) and satisfy 
^ol^Afi = h, and mq G if^(f2) be such that uo = F*t>o. Using formula 
with ip = u — Uq and -0 = f — fo, we see that m G Hg{Q) is a 
finite energy solution if and only if f = (^1,^2) £ H^{Mi) © H^{M2), 
u = F^v satisfies (EEl) for all G i^o (Mi) © H\M2), that is, satisfies 
equations flMl) and fl5Sl) in the weak sense. □ 

Assume next that A ^ spec(^i). Then the solution to equations ffTOj) . 
in the sense of definition fl36|) . may be found in terms of the resolvent 
7^l(A) of A, cf. ([32]). Indeed, comparing ([H]), ^ with equation fl55]) . 
we see that its solution u has the form 

(39) M = Uo + 7^l (A) (/ + • (Ti V^Mo - gwo + Amq) , 
wolan = h, supp(Mo) C ^] \ 5(2), 

at least when / G iv^(r^) and h G H^/'^{dQ) so that mq G i7^(fi). Since 

P(ai) = ©(^i^^) = i/o\g(fi), we see that Hg^{n) = V{A^^^^). There- 
fore, the operator 7^i(A) can be extended by continuity to a bounded 
operator from H~^{Q) onto Hq^IQ). This makes it possible to gener- 
alize (I39D for all uq G if](fi), i.e. h G H^/\dn), and / G //"H^)- 
Observe that the right-hand side in ( 1371) is related to the unbounded 
selfadjoint operator Aout in L'^i^), ^^(^oMt) C ifQ(n), associated with 
the form ai, see fll4l) where we use Q = Mi. More precisely, this 
operator is the unbounded selfadjoint operator in L^(fi) given by 

(40) Ao^t = -V^i.V'^^ + «:i, 



APPROXIMATE QUANTUM AND ACOUSTIC CLOAKING 15 

Moreover, when b G M'^), the selfajdoint operator associated to 

the form 0:2 on B{1) = M2 is the operator 

(41) = -Iv^^ ■ V^^ + «:2, 

V{A,n) = {V2 e H\B{1)) : d,v\9Bii) = 0}. 

When A ^ spec{Aout), the set (1571) coincides with the graph of the 
DN-map 

(42) ^out ■ '^ilan ^ {du + iJ^- (3i)vi\dn, 
where Vi solves equation (l34l) with /i = 0. 

Note that if 6 = 0, g = 0, then A = is an eigenvalue of ^1 with the 
corresponding eigenf unctions of the form (c/. ^3 ED]) 

^^^^ ""^""^-[co for a: €5(1), ^ 

However, as follows from Lemma l273l even, in this case the Cauchy data 
on dfl of solutions of the equation ( ITOl) with A = coincide with the 
Cauchy data of the solutions of Avi = on dMi = dQ. 

2.3. F— convergence and spectral convergence. In this section we 
establish F— convergence and spectral convergence results for as 
i? \ 1. To that end, following the comment in the previous section, 
observe that since a^-^ > a^^ for Ri > R2, one thus has 

(44) aR,[v,v]>ait,[v,v],veH\n), if > R2. 

This implies that R afi[v,v] is decreasing as i? \ 1. For R > 1, 
consider non-linear (quadratic) functionals : L'^{^) ^— > := M U 
{+00}, 

aR[v,v] when t> G //o(r2), 
00 otherwise. 



(45) aniv) 



For the ideal cloak, i.e., R= 1, define 

(46) a,(^) = l when G ffo\,(fi), 

1^ 00 otherwise. 

We will make extensive use of De Giorgi's F— convergence, see, e.g., 

mm- 

Definition 2.4. Let {Jr : 1 < R < 2} be a family of functionals 
on a Hilbert space H. We say that the Jr T-converge to Ji, or Ji = 
F-limjj^i J/j on U, if 
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(i) for every v & H, and all sequences vr converging to v in H as 
R\l, Jiiv) < liminfR^i Jr{vr); and 

(a) for every w & H there exists a sequence wr converging to w inTi 
such that Ji{w) > limsupjij^i JniwR). 

We then have the following result: 

Lemma 2.5. The functionals aR T-converge to ai as R\ 1, 

(47) T—\imaR = ai, on Ll{Q). 

Proof. By Lemma [2711 ifg^(f2) is a Hilbert space, when endowed with 
the norm (11^11^2 + '3i('u))^/^. Since the functionals Oo-^ are pointwise 

decreasing as i? \ 1, it follows from [161 Prop 5.7] that the functionals 
qr F— converge on to the 

(48) ai = sc^G, 

that is, the lower semicontinuous envelope of G (see [161, Def. 3.1]), 
where G : ^ M+ := M U {+00} is defined by 

(49) G(v) = I + qg^/M^) dx, for u G H^in), 

\ 00, otherwise. 

By [ini Prop 11.10], the function ai : Iv^(fi) — > MU {00} is a quadratic 
form. Moreover, by [HI Prop 12.16], its domain 

V{ai) = {ue L'l{VL) : ai{u) < 00}, 

endowed with the norm (||u||^ + 'ai{u)Y^'^ ^ is a Hilbert space. 

Now HliVL) is contained in both T'(ai) and T>{ai) = H^g^fl), and the 
norms of these Hilbert spaces coincide on Hq{Q). Moreover, by the 
proof of [201 Lem. 3.3], H^{n) is dense in V{ai). Thus P(ai) C V{ai). 

On the other hand, as G{v) > ai(f) for all v G L^{^) and ai is 

lower semicontinuous (see [IHl Prop. 2.16]), it follows that the lower 

semicontinuous envelope ai of G also satisfies ai{v) > ai{v). Hence 

V(ai) C V{ai). Thus, V(ai) = V{ai). □ 

Let us next consider the resolvent TZr{X) for A < 0. The following 
result then holds: 

Lemma 2.6. For any A < 0, the resolvents TZr{X), R > 1, strongly 
converge on L'^{^) to 7li{X), i.e., for any f G -^^^(^2), 

(50) lim7^^(A)/ = 7^l(A)/ 
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strongly in Lg{Q) and weakly in HQg{Q). 

Proof. The quadratic forms aji{u) — A||m||^, R> 1 in the Hilbert space 
Lg(f2) are associated to the unbounded selfadjoint operators Ar — \I. 
Thus we can use [Ml Thm. 13.6] to show that the resolvents {Ar — 
XI)~^ satisfy (|50ll in L^{fl). Indeed, to show the strong convergence 
f l50|) in L^{^) it is sufficient to prove the following three properties, 
(EI]), dSD and ([53]): 

(51) aR : Ll{Q) ^ M U {oo} are lower semicontinuous; 

(52) T—limaR = ai on Ll(Q); 

R\l " 

(53) ai{u) - \\\u\\l < liminf (aR^UR) - \\\ur\\1) as ur ^ u in Ll{fl), 

R\l 

where ^ denotes weak convergence in L^(f2). 

Clearly, the quadratic forms aR{u) — A||m||^ are lower semicontinuous 
on L^(n), see proof of Lemma [2751 By Lemma [231 the quadratic forms 
aR F— converge to oi. Thus ( l5Ti) and ( l52l) are valid. To see ( 153]) . it 
suffices to consider the case when ur^ u m. L"^ and 

lirninf (aR^UR) - \\\ur\\1) < oo. 

Next we consider a sequence \1 such that 

liminf (aR{uR) - \\\ur\\1) = lim (a\u'') - \\\u''\\l) , 

R\l ^ k^oo ^ 

where a^ = aRk and u'^ = URk G HliVL). Since converges weakly in 
Lg, the norms are uniformly bounded. Since also 

(54) ai{u^) < a'=(M^), 

we see that the sequence is uniformly bounded in Hg{Q). Let us 
now choose a subsequence of u'' (still denoted by u'^) which converges 
weakly to u in Hq^^Q). As the embedding ^ is compact, this 
sequence converges strongly in L'^{Q). Using the weak convergence in 
Hq^^Q), we see that 

(55) ai{u) - \\\u\\l = lim / fdiVV- + ^^/V''^ - A^^^Vm) dx. 
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By the Cauchy-Schwarz inequality, 



■ (^j {aiV\- V^u + g^'Wu\^ - WM") dx^ ' • 

Using (15^ and (!55|) . this imphes that 

aAu) - \\\u\\l < hminf (ai(u'') - A||m''||^)^^^ ■ (ai(u) - A||m||^)^''^ 

< hminf (akiu'') - A||m''||^)'^' ■ (ai{u) - A||m||^)'/' . 

The desired inequahty (153|) foUows immediately, proving (150|) in 

Finally, to prove the flSUl) holds weakly in Hlg{Vt), observe that, since 
<yR > 0-1, 

(56) / {crV^ ■ W + g^^'^q\v\'^) dx> I (aiV^ ■ + g^^'^q\v\'^) dx 
Jn Jn 

for V G Hq{Q). On the other hand, denoting ur = lZR{\)f and using 



A < 0, we have 



/ {(JrV^r ■ V'uR + g'^^q\uR\^) dx < 
Jn 



\ur\ 



The above two inequalities, together with the strong convergence (l50l) 
in L^(f2), show that the ||«_R||i/i are uniformly bounded. Thus, if weak 
convergence fl50|) in Hq^^Q) is not valid, there is a f 7^ m and a subse- 
quence Rn \ I such that, ^ w in H^. Then u^^^ — > f in L'^{^), 
which is in contradiction with the strong convergence (150|) in L'^{Q). 
Thus ([50D holds weakly in Hlg{Q). □ 

In some applications, e.g., dealing with scattering of plane waves e^<^'^^ , 
e M^, by the cloaking device, the main interest concerns not values 
A < 0, but rather A = > 0. To analyze this case, let us first 
consider the behavior of the spectra, spec(^ij), as i? \ 1. 

Lemma 2.7. Let K he a compact set with K fl spec{Ai) = 0. Then, 
for R sufficiently close to 1, K r\ spec {Ar) = 0. 

Proof. It suffices to consider K = [a, b] C M. Assume, to the contrary, 
the existence of a sequence Rn \ 1, /in ^ [O'yb], and functions G 
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Hi C Hlg, \\un\\g = 1, such that 

(57) AR„Un = /i„Mn- 

Then, 

7I/21 



\Un\ dx = fin < h. 



j ((^R^^^Un ■ V^Un + g^''^q\Un\^^ dx = fin J 9^ 

Therefore, as aR > ai and g > 0, this imphes llitnll h-i < h + 1. Thus 

rig 

there exists a subsequence of m„ and which we relabel as the original 
sequence, u G i^o ^ and ji G [a, 6], such that 

Un ^ u weakly in if^; Un ^ u strongly in L^; and /x„ — /i 

as n 00. Thus, in particular, \\u\\g = 1. 

Taking, e.g., A = —1 in Lemma [2.61 we know that TZr{—1) TZi{—l) 
as i? \ 1 in the strong operator topology on L^. Consider 7^i(— 
then 

7^l(-l)M = lim 7^K„(-1)M = lim 7^H„(-1)M„, in Ll{n) 

n— >oo n— >oo " 

where, in the last step we have used fact that A are non-negative 
operators yielding the estimate 1)||l2-*l2 < 1- 

However, taking into account (IFHl . 

7^/^„(-l)^^„ = — ^ — m„ m in Ll{Q) as n ^ 00. 

1 + /i„ 1 + /i ^ 

Thus, 7?.i( — !)« = (1 + fi)^^u with = 1, implying that /i G 

spec (.4.1) with u being an associated eigenfunction. This contradiction 
proves the statement. □ 

Lemma 2.8. Let K G C be compact subset such that Knspec{Ai) = 0- 
Then, for any f G L^{^) and X E K, 

(58) lim7^^,(A)/ = 7^(A)/ 

strongly in H^g^fl), and the convergence is uniform for X E K. 

Proof. Let 5 > satisfies 

Ks n spec(A) = 0, Ks = {zeC: dist(2, K) < 5}. 

It then follows from Lemma 12.71 that, for R sufficiently close to 1, 
Ks/2 n spec(^ij) = 0. As all Ar, R> 1 are self-adjoint in L'^{^), we 
see that, for R sufficiently close to 1, 

||^r(A)1|l^->l2 < ^, whenAGK 
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This implies that C \ spec(^i) = A;,, where Aj, denotes the region 
of houndedness for the family of operators Ar, R > 1, i.e., the set 
of those A G C for which the norms ||(^r — are bounded 

by some constant Ca > for all i? > 1, see [33i Sec. VIII. 1.1]. As 
spec(^i) is countable, A^ is connected. On the other hand, by Lemma 
12.61 R_ C As, where A^ is the region of strong convergence for the 
above family, i.e., the set of A G C such that for / G L'^{^) 



By [201 Lemma 3], H^{n) is dense in Hlg{n) =V{ai). Thus H^{n) is 
a core of the quadratic form ai[- , • ]. Now Hq{Q) = V{aji), for i? > 1 
and api are monotonically increasing with i? > 1 on HI{VL). Thus it 
follows from ^33* Th. VIII.3.6], that 



uniformly for A G -ft' where K is an arbitrary compact subset of A^. By 
Lemma [2.11 the desired convergence (l58l) now follows from (159|) - (!6T|) . 
□ 

Let IX ^ spec(^i). It follows from Lemma 12.71 that /i ^ spec(^ij) 
for i? > 1 sufficiently close to 1. For ^ spec(^i{), we denote by 
Nji{fj,) the subspace of L'^{^^) spanned by the eigenfunctions of Ar 
with eigenvalues A^ < /i. We also denote by P^, the orthoprojectors in 
Ll{Q) onto NRifi). By [221 Th. III.6.17], these (Riesz) projectors Pr 
have the representation 



where the contour F C C surrounds all the eigenvalues Aj of Ar satis- 
fying Xj < n and only those. 

Lemma 2.9. Let fi spec{Ai). For R sufficiently close to 1 



(59) lim7^^^(A)/ = 7^l(A)/ in L^n), 

it— >1 

Therefore, by [331 Thm. VIII.1.2], 

(60) A, = Ab = C\spec(A). 



(61) 



hmal(7^H(A)/-7^l(A)/) = 0, 

rt\l 



(62) 




(63) 



dim{NR{^)) = dim{Ni{fi)). 



Moreover, 



(64) 
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Proof. Recall that 'D(afi) = Hq^VI) are independent of i? > 1, -f^o (^^) C 
T'(ai) = Hq^^Q), and afi[u,u] are decreasing, as \ 1, for all u G 
vlaR). Thus the identity §^ follows directly by [33, Th. VIII.3.15]. 



Using representation fl62|) and Lemma [2.81 we see that 
hm = Pf strongly in L'^{i^)- 

As and Pj" are orthoprojectors, this and (l63ll yields (Ell) by |33] (see 
Lemmas VIII.1.23 and VIII.1.24). □ 

We remark that in the course of this paper we need a number of results 
concerning convergence of orthoprojectors that appear similar to ( IMj) . 
e.g., ( !78|) . (11141) and (11331) . but these are for different operators or with 
respect to different operator norms and require separate proofs. 

2.4. Approximating the singular bulk modulus g^^'^ by nonsin- 
gular densities. Above, in the operator Ar = —g~^/'^V^- orV^ + q 
there appears the determinant g of the metric (also denoted (7!), which 
vanishes at the cloaking surface S. We now consider how to approx- 
imate the scalar function g by functions gm that are bounded from 
below with positive constants. To this end, we introduce the functions 

(65) (7m(a;) = max (5f(x), 1/m) , m G Z+. 
Then L\n,gll^dx) C L2(fi) and 

(66) 11/11, < 11/11^.(^,^1/2,,), for / G L'in,gl/'dx). 

The multiplication map / ^-^ g^^'^f is unitary from = L'^{^1, g^dx) 
onto L'^{n,g-^dx). Note that L^{n,g--^dx) C /.^(fi) c Lg(fi). Next 
we will consider operators g^^'^ARU. For / G L^(f2), we have 

(67) {Ar - X)u = /, 

where both sides are in L^(f2), if and only if m is a solution to the 
acoustic equation with mass density tensor aR, bulk modulus g^^"^, and 
potential qg^^'^), 



(68) {g'/'AR - Xg' 



/2^ 



u 



where F = g^^'^f G L'^{Q,g ^dx). By the above considerations, we 
have for / G L1{VL) 

(69) {Ar - A)-V = {g'^'AR - W^r^g'/'f). 
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Later in this section we keep R > 1 fixed. Define an unbounded self- 
adjoint operator Br in L^(f2), having the same differential expression 
as the operator q^^'^Ar, but with different domain, 

(70) Bnu := -V" ■ (anV'u) + g^'\u, 

V{Er) = {ue Hl{n) : • {(JrV'u) e L\Q)}. 

Since V{Ar) C Hq{Q), see (I5U]) . Br is an extension of q^^'^Ar and, 
in particular, ■ (ctrV^m) G L'^{Q, g~2dx) for u G V{Ar); however, 
■ (a/jV^n) G L2(fi) for m G V{Br). Note that by ([69]), 

(71) iAR-\r'f = iBR-W/'r\g'/'f) for/GLJ(^]), 

where A ^ spec(^/j). We will use this formula extensively later for 

/ G L^{n) c mn). 

Lemma 2.10. The operator Br — Xg^^"^ has a bounded inverse if and 
only if \ ^ spec (Ar) . 

Proof. For A < the operator Br — Xg^^'^ is strictly positive and, 
since T>{Br) C HI{VL), has a compact resolvent. Therefore, the oper- 
ator [Br — Xg^^"^)^^ exists for A < and is bounded in L^(f2). Since 
the multiplication, u i— > g^^'^u, is bounded in L^(f2), by the analytic 
Fredholm theory [56] the operator-valued function 

X^{Br- Xg'/')-' = [/ - (A + 1){Br + g"')-'g'/'] ' {Br + g^l^Y' 

is a meromorphic operator-valued function of A G C. Therefore, if the 
inverse {Br — Xg^^"^)'^ does not exist for a given A > 0, then there is 
V G HI{VL) such that 

{Br - Xg^'^)v = 0. 

In this case Brv = Xg^^'^v G L'^{^1, g~^dx) and we see that v G ^^{Ar), 
i.e., X G spec(^/j). On the other hand, if A G spec(^ij), i.e. 

-g-^/^V^ . {aRV^u) +qu = Xu, 

then V^- (afiV^'M) = qg^/'^u - Xg^'^u G L^(fi), z.e., m G V{Br), so that 
Br — Xg^^'^ does not have a bounded inverse. □ 

Next we consider the uniform convergence of resolvents. To this end 
we introduce operators Br^^^ m G Z+, in L^(fi), of the form 

(72) BR^mU := -V" ■ {(JrV^u) + qg]l\, 

V{Br,„,) = {ue H',{Q) : V' ■ {aRV'u) G L\Q)} = V{Br). 
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The operator BR^m is associated with the operator Ar^ju in the same 
way that Ar is associated with Br, where the operator AR^m is the 

self-adjoint operator in L'^{Q, gU^dx) defined by 

(73) AR,mU := -(7~'/'V' ■ {(TrV'u) + qu, 

V{AR,m) = {ue H'oin) : ■ {aRV'u) E L\n)}. 
Note that V{Ar,J = V{Br^^). 

Lemma 2.11. Let R > 1 and K G C be compact and such that K fl 
spec (Ar) = 0. Then there is an rriR G Z+ such that KHspec (AR^m) = 
for m > mR, and 

(74) hm II {Ar,^ -X)-^- {Ar - ||L2(n)^//i(n) = 0, 
uniformly for X E K. 

Proof. By the assumptions on K, the inverse {Br — \g^^'^)~^ exists and 
is a continuous function of A G -ft' with respect to the L^(f2)-operator 
norm topology. Let 

d = max\\{BR - \g^^'^)~'^\\L^(n)^L2{n) < oo. 

Denote VR^m{\) = (A - q){gll^ - g^''^) G L°°(^^), so that 

(75) lim ||V}?,„(A)|Uoc(f,) = 

m— >oo 

uniformly for \ E K. Thus there is m{K) > such that || VR_m(A)||L°° < 
{2d)~^ for m > m{K), X E K. Therefore, {BR^m — Af^m^)"^ exists for 
X E K and is given by 

(76) {Br,^ - Xg'fy' = {Br - Xg'/y' [l + Vr,U^){Br - Xg^/^T 

where the right hand side can be written as a Neumann series. This 
also shows that there is an niR such that K fl spec {AR^m) = for 
m > niR. 

For any A' G K, {Br — X'g^^"^)'^ is a bounded operator from L^(fi) to 
Hq{Q), and for if |A — /i| < {2d)~^ we have 

{Br - Xg'/y' = {Br - X'g'/y' [l + (A' - X)g'/\BR - X'g'/Y'] ' ■ 

Using this we see that the norm of {Br — Xg^/'^Y^ : L'^{Vl) — > Hl{Vl) 
is uniformly bounded m. X E K. Using formulae (175|) and (!76|) . we see 
that 

(77) lim II {Br,^ - XgllY' - {Br - Xg'/')'' \\LHn)^H^,in) = 
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uniformly for A G -fT, and that the norms of operators {BR^rn — ^gm^) ^ : 
L^(n) Hq{Q) are uniformly bounded for X & K. This proves 

lim II {Br^^ - XglC^) ^ - {Br- Xg^^"^) ^ |L2(n)^m(f^) = 0- 

Additionally, the multiplication operators gW^, g^^"^ are bounded on 
L^(f2), uniformly in m, and 

llfl'm^ - g^^'^Wh^^L^ < m'^^'^ ^0, as m oo. 



Together with equation fl77j) and the boundedness of {BR^m — ^g 
as operators from L'^{Q) to Hq{Q), this implies that 



1/2. 
m I 



lim 

m— >oo 



{Bn,m-XgUT'9U'-{BR-Xg'/Y'g 



1/2 



= 0. 



This in turn implies equation fl74p . due to formula fl71l) and the relations 
between BR^m, Br and AR^m, -Ar, which follow from their definitions 
dsn]), dTUHZi, and ([73]). ' □ 

Let R > 1, fi ^ spec(^ij,). It follows from Lemma 12.111 that fi ^ 
spec(.4i?,rn) for m sufficiently large. Denote by NR^mifJ^) the subspace of 
L'^{Q, gW^dx) spanned by the eigenf unctions of AR,m with eigenvalues 
in (— oo,/i), cf. the definition of NR^fx). Also denote by -P^^ the 

orthogonal eigenprojectors onto NR^rn{lA in L'^{^, gW^dx). 

Clearly, since, for m G Z+, L'^{Q,, gli'^dx) = L^(fi) as sets, we can con- 
sider Pr^ as projectors, although not orthogonal, on L^(f2). Recall 
that is an orthoprojector in L'^{^) onto NR{fi) C Hq{Q). Restrict- 
ing it to L^(f2), we obtain a projector, which we still call Pr, on L^(fi). 
Again, P^ is not an orthoprojector on L'^{Q). However, we can compare 
these projectors, as well as spaces NR^m^lj), Nr^h). 

Corollary 2.12. Let R > 1 and /i ^ spec{AR). The Riesz projectors 
^Rm corresponding to the operators AR^m satisfy 

(78) lim ||Pi^,„-Pi^|U.(^)^^i(^) = 0. 
Moreover, for m sufficiently large, 

(79) dim{NR^rn{li)) = dim{NR{fi)). 



Proof. Let F C C be a contour surrounding only Aq from spec(^/{). 
When m is large enough, the Riesz projectors -Pr^ have a representa- 
tion analogous to (l62l) . obtained by replacing Ar by AR^m and using 
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the contour T. Thus 



hm ll-P^m - -PR||L2(f7)-*Hi(n) < 



TTi — >c>o p 



where dl is the arclength measure on T. Taking into account ( 174|) . this 
formula imphes ( 1781) . 



Using equation (1781) . we see that there exists an mo such that for m > 
mo we have ||P^,„ - PR||L2(Q)_,L2(f^) < 1. Using, e.g., [33l Cor. IV.2.6], 
we see that this proves (179|) . □ 



3. Approximating anisotropic by isotropic mass densities 



We now show, using techniques from homogenization theory, cf. IH 
[131 [16], that we can approximate arbitrarily closely, on the level of the 
operators, the nonsingular anisotropic approximate mass densities a^, 
for any fixed i? > 1, by a family of nonsingular isotropic mass densi- 
ties cTi?_e, e > 0, which will thus also function as approximate cloaks. 
This can be considered as the reverse of the traditional homogenization 
theory. 



3.1. Inverse homogenization with magnetic potential. Observe 
that, as all the approximate cloaks ctr, R > 1, are rotationally invari- 
ant, it is natural to use spherical coordinates. Namely, we will use either 
the Euclidian coordinates x = {x^, x^, x^), or the spherical coordinates 
s = s{x) = {r{x),9{x),(f){x)). Note that we use the same notation, 
X, for a point inside Q and its Euclidian coordinates, x = {x^,x'^,x^). 
Which meaning is intended will be always clear from the context. We 
denote by X : {x^,x^,x^) i— > {r,9,(f)) the corresponding coordinate 
transformation. To exploit the rotational invariance, we will employ in 
the homogenization process, cf. [H [13] , cells which are cubes in these 
spherical coordinates. 

To approximate the anisotropic mass densities in spherical coordinates 
in the ball fl = 5(0, 3), let us consider isotropic mass densities of the 
form 

(80) (Xeix) = cr{x,r{x)/e), a{x,r') = h{x,r')I e R^''^ 

Here h{x,r') is a scalar valued function, to be chosen later, that is 
periodic in r' with period 1 and is bounded from above and below, i.e., 

(81) h{x,r' + 1) = h{x,r'), < Ci < h{x,r') < 02- 
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We will consider densities for which h{x, r') is independent of r' for x 
with r{x) < 1 and 5/2 < r(x) < 3, that is, 

(82) h{x,r') = h{x) if |x| < 1 or 5/2 < |x| < 3. 

We make this assumption since later we will use the isotropic mass 
densities to approximate the non-singular anisotropic mass densities 
ctr, R> 1 that are isotropic for r{x) < 1 and 5/2 < r(x) < 3. 

Let (r, 9, (p) and (r', 6' , 0') be spherical coordinates corresponding to 
two different scales. Then, in these coordinates, 

(X,(a,))(s)=a*(s,s/£), 

where 

a*{s,r') = h{s,r')dei{DX{x))-'DX{x)DX\x)\x(.)=s 

I r^sin^ \ 
= h{s,r') sin^ 

I / 

Here and later we denote by a, sometimes with various indices, various 
mass tensors (or matrices) in the Euclidian coordinates, while a* always 
stands for their representation in the spherical coordinates. 

In the following, the material on homogenization is a quite straight- 
forward generalization of known results [H [13] . However, as we need 
to introduce changes due both to the presence of a magnetic potential 
and the use of spherical coordinates, for completeness we give details 
of the arguments. 

In the small-scale coordinates t = (r', 9'), we denote by = (1, 0, 0), 
= (0, 1, 0), and = (0, 0, 1) the vectors corresponding to the differ- 
ential forms dr', dO' and dcj)' , respectively. Let W^{s,t), j = 1,2,3, be 
the solutions of 

(83) Vfcr*{s,t)i\/tW\s,t) + e') = 0, t = (r', 0') ^ M^ 

that are 1-periodic functions in r', 9' and 0' variables (noting that the 
periodicity in r', 9' and 0' has no relation to periodicity in the "large- 
scale" spherical coordinates 6', 0), and satisfy, /^q W\s^t)dt = 0, for 
all s, where dt = dr'd9'd(f)'. Since a* is independent of 9', 0', the above 
conditions imply that = for j = 2, 3. As for VT^, it satisfies 

d_ f _ dh{s,r') 
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with being 1-periodic with respect to {$', 0'). These imply that 
is independent of {9', cj)') with 

dW^ ^ Co 
-1 + 



dr' h{s,r') 

To find the constant Cq we use the periodicity of W^, now with respect 
to r', to get that Cq is given by the harmonic means h'^"-^"^ of h, 

dr' 



U) Co = : 



1 j^' ^ 



Ks,r') 



Define the corrector matrices [Ij as 

(85) P^{s,t) = -^W%s,t) + 5^. 

Then the homogenized mass density in the spherical coordinates, cr^o^, 
is given by 

(86) ialJ^' (^) = E / ^K'(-'^ ^) dt. 

Z~1 i [0,113 



p=i 



We note that, applying integration by parts and using definition 
and equations (!83l) . formula (l86l) can be written also in a more sym 
metric form 



p,q=l 



Using the above formulae for the W\ it follows from fl86|) that 
/ h'"'''"'{s)r^sm{e) 

\ /i"(s)sin-^(0) 

Here h^-is) denotes the arithmetic means of h in the second variable, 



h^{s) = / h{s,r')dr'. 

J [0,1] 

Returning to the Euclidian coordinates, one sees that the conductivity, 
(^hom{x) = X^al^^{x), has the form 

(87) CThomix) = UJi{x)n{x) +UJ2{x){I -n{x)), 

with 

(88) tu^ix) = /i'*""'"(x), uj2ix) = h\x) 
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and n(x) : M'^ — > M'^ being the projection to the radial direction, 



Il(x) V = { V ■ 



x\ / \x\ 



represented by the matrix (|x| '^x^x'')^ i^^-^. 

Next, we analyze the Dirichlet problems for elliptic equations 

(89) - V' ■ fT, V V + Que = /, u,\dn = h. 

Here, h{x) = 62(3^)) ^3(3^)) is the magnetic potential and Q{x) 

is a scalar function, with Q,bj G L°°(f2;]R), and as are isotropic mass 
densities bounded from above and below by positive constants indepen- 
dent of e. Moreover, / G H~^{il), and h G H^^'^{dVL). Our goal is to 
show that the solutions Ue convergence to the solution of the equation 

(90) -V^ ■ahomV^u + Qu = f, u\an = h. 

By adapting the technique of Allaire pj, we can prove the following 
result. 

Proposition 3.1. Let (Te, e > he mass densities in f2 satisfying ^8D^, 
[8l\). and [8^) . ahom be the mass density defined by ( [g7| j, Q G L°°(f2), 
Q{x) > 0, andb= (61 (x), 62(2;), 63(3;)) be a vector field, b G L°°{Q;R^). 
Then the solutions Ue of equations and solution u of equation ^9D\) 
satisfy 

(91) limw^ = u weakly in H^iVt). 

Proof. Let Eh G H^{^), supp{Eh) C {5/2 < \x\ < 3} be an extension 
of h, i.e., Eh\dn = h. Writing = E^, + and u = E^, + v, we see 
that the functions Ve and v satisfy 

(92) -V' ■ aeV\ + Qv, = f, on Q, v.lan = 0, 

(93) - ■ cXhomV^ + Q^; = / onfi, vlon = 0- 
Here 

(94) f = f-^K a.v'Eh + QE^ G H-\n), 

is independent of e. For proving it is enough to show that 
converges to v weakly in Hq{Q). 

Let us recall that by Lemma [Z2| is (strictly) positive definite. Since 
Q{x) > 0, we see that there exists a cq > such that 

{\V\{x)\'^ + Q{x)\u{x)\'^) dx > co||M||^i(f^) for u G H^{n). 
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As the mass densities cxe are uniformly bounded from above and below 
by positive constants, it follows from the proof of Lemma [22] that there 
is also Ci > such that 

(v^- aeVhi + Q{x)\u\'^^ dx > Ci\\u\\]^^^^^ for u e H^{^). 

Thus, using the Lax-Milgram lemma we see that the solutions ^^(x) of 
equations (l92l) satisfy 

(95) \\ve\\Hl,{n) < c\\f\\H~Hn), 

where c > is independent of £ > 0. Therefore, the solutions Vs{x) are 
uniformly bounded in Hq{Q). Thus, for an arbitrary sequence —>■ 0, 
the corresponding v^^ have a subsequence that converges weakly to 
some function w{x) in Hq{Q). Let us show that w coincides with the 
solution V of equation flU^ . 

To this end we consider convergence in a finer, two-scale, sense on local 
coordinate neighborhoods. Let U C be an open set on which we 
can define, in a regular manner, spherical coordinates. For example, 
by choosing two antipodal points as the South and North poles and 
connecting those by a meridian 7, we can take U so that ?7 CC §^ \ 7 
and define polar coordinates on U. Let Q' = {ru : G f/, r G 
{fi,f2)} C Q with some r2 > ri > 0. Clearly, X : fi' ^ defines 
the spherical coordinates, x ^ {r,6,(j)), with domain W := X{Q') = 
(ri,r2) X U. 

Rewrite now the equations (l92l) and (l93l) on W in these spherical coor- 
dinates and multiply the equations so obtained by det{DX{X~^{s))). 
The resulting equations are 

(96) -(V, + ib*{s)) ■ a*{s, '-){Vs + ^b*{s))v,{s) + Q*{s)ve{s) = f{s), 

(97) -(V. + tb*is)) ■ al^is)iVs + tb*is))vis) + Q*is)vis) = /(s), 
where s & W. Here 

Q*{s) = det{DX{X-\s)))Q{X-\s)) e L^{W), 
b*is) = d,X\X'\s))bkiX~\s))eL'^iW;R''), 

correspond to the electric potential Q and magnetic potential b in the 
spherical coordinates 

and f{s) = det{DX{X-\s)))f{X-\s)) G H-^{W). To simplify no- 
tations, we continue to denote the functions v and Vg in the spherical 
coordinates by v{s) = {{X~^)*v){s), Ve{s) = {{X~'^)*Ve){s). Finally, 
Vs is the vector field {dr, do, d^). 



/ 

Jn 
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In the following, even though in the forthcoming applications the mass 
densities cr*{x, r{x)/e) will depend only on the small-scale variable r' = 
r{x)/e, we consider the general case when the mass densities depend on 
all small-scale coordinates t = (r',^',0'). Let T = [0, 1]^ and C^{T) 
denote those C^iT) functions that can be continued as Z^-periodic 
functions in which are in C"^{E?). By definition, a family 
Ve G L'^iyV) is said to two-scale converge, as e 0, to a function 
Wo(s, t) G L2(1^ X T) if, for all test functions i){s, t) in C^{W] C^{T)), 
we have 

(98) lim / Ve{s)ilj{s, -) ds = / / Vo{s,t)ilj{s,t) dsdt. 

Jw ^ Jw Jt 

By [1], the two-scale convergence of Ve implies the weak convergence 
of Ve{s) in L'^iW) to the function w{s) = JrpVo{s,t) dt, so that the 
two-scale convergence gives finer information on the convergence than 
the weak convergence. For example, functions of the form u{s,s/e) 
two-scale converge to u{s,t). 

By [I], every bounded family ^^(s) G L'^{W) contains a two-scale con- 
verging sequence. Moreover, if fe(s) is a bounded family in H^{W) 
that converges in L'^{W) to w{s) as e — » 0, then Ve{s) also two-scale 
converges to w{s) and there is a function wi{s, t) G L'^iW x T), so that 
Vsfe(s) two-scale converge to V sw{s) + VtWi{s,t). For example, if 
functions Ve{s) have the form Ve{s) = uo{s) +eui{s, s/e), then Vsfe(s) 
two-scale converge to u{s,t) = V'Uo(s) + VtUi{s,t). 

As noted above, the solutions ^^(x) of equations fl92l) are uniformly 
bounded in Hq{Q). Consider the restrictions of these functions on Q' 
and rewrite them in the spherical coordinates on W. Then any sequence 
of fe(s) has a subsequence Ve {s), j G Z^, Sj ^ as j ^ oo, that 
weakly converges, in H^{W), to some function w{s). By ^ v^.{s) also 
two-scale converge to w{s), as e — *• 0, and there is a function Wi(s, t) G 
L'^iW X T) so that Vv^^^s) two-scale converge to Vsif(s) + V tWi{s,t). 

Let us now multiply both sides of equation flUBl) by a test function 
0(s)+£0i(s, s/e), where 0(s) G Co~(W^) and 0i(s, t) G {W ; {T)) , 
and integrate over VT. Using integration by parts, we obtain 

(99) / a*(s, -)[(V. + ih*{s))vM] ■ iViis, -) + eV2{s, -)) ds + 

Jw ^ 6 6 

+ [ Q*{s)v,{sms)+6Ms,-)]ds= [ f{sms) + 6Ms,-)]ds, 
Jw ^ Jn' ^ 
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where 

r]i{s, t) = V.0(s) + ih*{s)(l){s) + Vt0i(s, t), 
m{s,t) = VsMs,t) + ib*{s)(f)i{s,t). 



Substitute e = ej in ( l99l) and let j ^ oo. First, we observe that since 
ri2{s, s/e), e G (0, 1) are uniformly bounded in L'^{W), the terms in fl99|) . 
involving 172, tend to zero. Second, as Ve{s) converge in strong topology 
of L'^{W) to w, the first integral in fllOOp tends to J,^ Q*{s)w{s)(f){s) ds. 

Third, the inner product of f{s) and e(j)i in the last integral in fllOOl) 
goes to zero. Fourth, the functions ipsi^s) = a*{s,s/e)rii{s,s/e) two- 
scale converge to the function il^{s,t) = a*{s,t)rii{s,t). Since ip^x.y) G 
L2(fi;C^(r)), we have by [H Lem. 1.3] that 

(100) lim||?/'£||i2(H/) = ||^||l2(h/xt)- 

Furthermore, as Vfe.(s) two-scale converges to V sw{s) + V tWi{s,t), it 
follows from ffTOOD and [H Thm. 1.8], that 

/ a*{sA)[{ys + ih\s))ve^{s)]-n,{sA )ds^ 

/ / a*{s, t)[Vsw{s) + VtWiis, t)] ■ r7i(s, t) dsdt, 
Jw Jt 

as j 00. Summarizing, we see that 

(101) / / a*is,t)[Vsw{s) + \/tWi{s,t)+tb*{s)wis)] 
Jw Jt 

■ [V,0(s) + ib*{s)(p{x) + (Vt0i)(s, t)] dsdt + 



+ / / Q*is)w{s)(l){s)dsdt= / f{s)(j){s)ds. 
Jw Jt Jw 

Taking (j){s) = in fllOll) and varying 0i(s,t) over all test functions, 
integration by parts with respect to t shows that 

(102) Vt ■ a*{s,t)[Vsw{s)+VtWi{s,t)+ib*{s)w{s)]=0 a.e. in x T. 

Since V sw{s)+ib*{s)w{s) depends only on s and thus can be considered 
as a constant vector in the t-variable, we see that 

(103) wi{s,t) = J2 (wjis)+ib*is)w{s)] W\s,t) a.e. inW xT, 
where are solutions of the cell equations (!83l) . 
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On the other hand, taking 0i(s, t) = in fllOll) and varying 0(s) over 
all test functions, we see, using integration by parts with respect to s, 
that 

-{Vs + ib*{s)) ■ (^j a*{s,t)[V,w{s) +VtWi{s,t) +ib*{s)w{s)]d?j + 

(104) +Q*{s)w{s) = f{s) a.e. in W. 

Comparing the above integral with formula (186!) for ahom, we see that 
( 11031) and (ITOD imply that 

-(V. + tb*is)) ■ <_(s)(V. + tb*{s))w{s) + Q*is)w{s) = f{s) in W, 
or equivalently, in the Cartesian coordinates 

(105) - (V + ib) ■ ahomiy + ib)w + Qw = f mQ'. 

Since Q \ {0} can be covered with coordinate neighborhoods Q' used 
above, the equation ( 11051) is valid in the domain Q \ {O}. Moreover, 
since all the mass densities ahom and are the same near the origin we 
see that all functions v^^ (x) satisfy the equation (11051) near the origin. 
As Vej{x) converge weakly in -f^o(^^) to w{x), the equation ( I105p is 
valid in Q. This means that w{x), as well as v{x), is the solution of 
the equation (jHSD- By the Lax-Milgram Theorem, equation (1^^ has a 
unique solution, and thus w = v. 

Summarizing, we have now shown that an arbitrary sequence of the 
original family Vs{x) has a subsequence that weakly converges, in Hq{Q), 
to the solution v{x) of equations (193|) . This shows that the whole fam- 
ily Ve{x) of the solutions of the equations (192|) converge weakly to v{x) 
in H^{n). □ 

3.2. Approximation of the mass densities ct/j. Next we apply the 
above results to approximate the mass densities aji for a fixed R > 1. 
In the forthcoming analysis, we will choose the function h to be 

106 h(s,r') = , ; / , , , 

where p{r') is a fixed positive, smooth 1-periodic function of r' and 
a{s),b{s) are positive. For any x & Q and a given function p{r'), the 
pair {uji{x) , u!2{x)) in ( 187|) . ( IHHl) depends only on the values of functions 
{a{s{x)),b{s{x))). By choosing (a(s),6(s)) appropriately, it is possible 
to achieve any value 

{tOi{x) , uj2{x)) G M^, with tOi{x) < uj2ix). 
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We can choose a{x) = aR^x) and h{x) = hR^x) so that o-hom{x) corre- 
sponds to the approximate cloaking mass density, that is, 

(107) (yhom{x) = o-r(x), X en. 

In sequel, we will denote by aR{x,r') = hR{x,r')I the isotropic mass 
densities for which the homogenized mass densities ahom satisfy fll07p . 
Observe that, for i? > 1, we can choose hR{s,r') so that 

(108) < Ci{R) <hR{s,r') <C2, Ci{R) ^ as R\ 1. 
Let A < and BR^m,e be the operators 

V{BR,„,,e) = {ue H',{n) : -V'-aR,eV'u G L\Q)}. 

Lemma 3.2. For any A < and f E L'^i^), 

\im{BR,m,s - XgllY'f = i^R^n - Xgll'r'f, 

where BRm is as in (7^ ), and the limit holds both in the sense of the 
strong topology in L^(fi) and weak topology in H^{Q). 

Proof. This follows immediately from Prop. 13.11 by taking Q = {q — 
X)gll\ □ 

As in the previous section, one can analyze the convergence of the 
resolvents in more detail. To this end, introduce operators AR,m,e = 

gm^'^BR^rn,e that is, 

(109) AR,m,e = -^7™ /'V'' • (T«,,V^n + qu, 

'DiAR,m,e) = V{BR,^,e) = {u E HH^I) : -V'- (T^,,V'm G L2(fi)}. 
Clearly, the operators AR^rn,e are self-adjoint in L^(i7, gW^dx) and 

(110) {Ar,^,, - x)-'f = {BR,m,e - ^gli^r\gll'f) 

for A ^ spec {AR,m,e) and f e L^Q), cf. m, 

Lemma 3.3. Let R > 1, m E Z+ and K G C be compact and such 
that dist{K, spec{AR^m)) > <5o > 0. Then 

(i) There is an = eo{K, R, m) such that such that for < e < Eq 

2 

\\{AR,m,e - ^y'\\L^^n,gll'dx)^LHn,gU'd.) ^7"' ^ G /T, 



(a) For any f G L^{^), 

(111) lim(^K,^,, - A)-V = {Ar,^ - A)-V 
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in the strong topology of L'^iVt) and weak topology of Hq{Q), with con- 
vergence being uniform for X E K. 

Proof, (i) Since the AR^rn,e are self-adjoint operators on L^(f2, gW^dx), 
it is enough to show that there is £o > such that the operators 
•^R,m,e — are invertible for |A — /i| < (5o/2 if £ < ^o- 

Assume that there are \ 0, /i„ — > fi, |A — yu| < So/2, and u„ , with 

uch that 



\''-\\LHn,gl/'d.) = 1' such that 



Rewrite this as 

(112) - V- CTR^^^VUn = fn ■= -gli'^QUn + /inS'm^Mn, Un\dn = 0. 

Taking into account the uniform boundedness of the right hand side of 
ffml) in L'^{n), it follows from ([92]), ([95]), with Q = 0, that < 
C, for some C > independent e„. 

Thus, up to a subsequence, there exists u G Hq{Q) such that 

n — >■ cx) in the weak topology of H^{Q) and strong topology of L'^{Q). 

This shows, in particular, that il'^llj;^2(f^ j,i/2^^-) = 1- 

Let us show that u is an eigenfunction of AR^m, see ( 1731) . corresponding 
to the eigenvalue /i; as |A— /i| < 6o/2 with \ E K and dist(i^, spec(^i?,m)) > 
6 > 0, this would yield a contradiction. We compare Un with Vn, the 
solution to 

(113) - V^- (TR,e„ W = 7„ := -gli\u + figl/\ Vn\an = 0. 

Letting /i„ — /i, (11131) and the weak convergence m„ — > m in iJ^(fi) 
imply that 

lim II /„ - fn\\L2{n) = 0. 

n— >oo 

Appealing again to (llOSp . with R > 1 fixed, we see that m„ — v„ — »• 
in the strong topology of Hq{Q). 

Using Proposition 13. II with Q = 0, we see that f „ ^ f in H^{Q), where 

Summarizing, we have that u = v, showing that u is an eigenfunction 
of Ar corresponding to the eigenvalue A. This proves claim (i). 

(ii) Compare the solutions of 

-V^- CT/{,eV^Mm,£ = fm,e ■= f - 9m'^qUm,s + \g]i^U^^e, Ura,e\dQ. = 0. 
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with those of 

-V^- anV'u = f^:=f- gK^u + Xg'J^, uU = 0. 

It follows from claim (i) together with fl95l) . that the Um,e are uniformly 
bounded in Hq{Q). Choosing a weakly converging subsequence, we 
show, by similar reasoning to the above, that the limit is Um- This 
proves 

e— »0 

This, f ill 01) . and the boundedness of the multiplication operator / i— > 
gH^f in L\n) yields dm]). □ 

For /i ^ spec {AR^m,e) we denote by Nji^rn,e{f^) the subspace of L^(fi) 
spanned by the eigenfunctions of AR^m,e corresponding to the eigen- 
values in (— oo,/i), and by Ppime the orthoprojectors in L'^{Q, gll'^dx) 
onto A^i?„m,e(/i). 

Lemma 3.4. Let R > 1, m E and ^ ^ spec{An^rn)- Then 

(i) There is an SR^m > such that /i ^ spec {AR^m,e) for < e < SR^m- 

(a) For f G L^(fi) we have 

(114) lim Pi^,^, J = Pi^,^/, 

in the strong topology of L'^{Q) and weak topology of Hq{Q). Moreover, 

(115) hmdim{NR^rn,sil^)) = dim(A^R,„(/i)). 

Remark. We note that the strong convergence of the resolvents does 
not generally imply the stability of the resolvent set, cf. [33l Sect. 
VIll.1.2]. 

Proof, (i) Take K = {fi} and Sq = dist(ii', spec(^R_m))- Lemma [3731 

(i) implies that, for sufficiently small e, dist{K,spec{AR^rn,e)) > ^o/2. 

(ii) Relation flll4p follows immediately from fill II) and the Riesz for- 
mula for the projectors, cf. flB21) . For the proof of 01151] . let us 
consider Aq G spec (AR^m) and a,b E R, a < Xq < b such that 
[a,b] n spec (AR^m) = {Xq}- Then by (i) there are 0{e), 9{e) 
as e — > such that 

spec iAR,m,e) n [a + e{e), Aq - e{e)) = 0, 
spec iAR,m,e) n (Ao + 9(6), b - 9(6)] = 0. 
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Let = P^;^y - P^+„'J) and = P^^ - P^^. Then, by (HH, 
linig^o-P'^/ = in the strong topology of L^(f2) and weak topology 

To establish (11151) . it is sufficient to show that 

(116) limdim(Ran(P")) = dim(Ran (P°)). 

e— >0 

Indeed, if we show (11161) for each eigenvalue of An^m smaller than /i, 
equality (IllSp follows. 

To show (I116p . we note that the ranges of the orthoprojectors P*^, P° 
are also subspaces of Hq{Q) C L'^{Q). We show first that 

(117) liminf dim(Ran(P^)) > dim(Ran (P°)). 

On the contrary, if this does not hold, there is a sequence e{n) — > 
such that 

(118) dim(Ran (P"^"))) < dim(Ran (P°)) = ko. 

Denote by 77^, k = 1, . . . ,kQ, an L^(r2)-orthonormal basis in Ran (P*^). 
(Note that as Ran (P°) C Hq{Q) is finite dimensional all norms are 
equivalent). Introduce rjn^k = P'^Vk ^ Ran(P^), k = l,...,kQ. Con- 
sider the Gram-Schmidt matrices, G" = 



Gki = / Vn,k{^)Vn,iix) dx. 
Jn 

By (11141) and orthonormality of rjk, the matrix is invertible for 
sufficiently large n. Thus, rjn^k, k = 1, . . . , fco are linearly independent. 
This implies that dim(Ran (P^*^"-*)) > ko, contradicting (11181) . This 
proves (I117p . 

Assume next that 

(119) limsup dim(Ran (P^)) > dim(Ran (P°)). 

Then there is a sequence e{n), such that 

dim(Ran(P"("))) > dim((Ran (P°)). 
Thus, there are such that 

(120) n„GRan(P^(")), = I, 

(121) (m„, v)^2^^^gi/2^^^ = 0, for all v G Ran (P°). 
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Then 

Mn = ^ Un,ki^n,k, ^ \un,k\^ = 1, K = dim(Ran (P^^"^)), 

k=l k=l 

where ijjn,k are 5fm^c?x)-orthonormal eigenfunctions of AR^rn,£{n), 

•y^R,m,£{n)tpn,k = ^n,ktpn,k, \^n,k — < d{e{n)) 

Therefore, 

kn 

(122) AR^rn,e{n)Un " AoMn = /„ := ^Mn,fc(A„,fc - \o)^n,k ^ aS W ^ OO 

k=l 

in L'^{VL, gll"^ dx), so that 

= (/n, «n)^2(f,,3V2rf,) ^ as n ^ oo. 

Together with fll20p and the fact that q G L°°{VL), this imphes that the 
forms crij,£(n)V*Mn ■ V^Undx are uniformly bounded in n. Recalhng 
f l95|) . we see that with some c > 0, 

< c, for all n E Z+. 

Restricting to a subsequence of the assume there exists uq G 

H^{fl) such that 

(123) lim Un = Uq weakly in Hq{Q) and strongly in //^(fi), 

n— ►oo 

SO that also ||^o||j^2(q gi/z^^-, = 1- Let us show that An^mUo = ^qUq,, con- 
tradicting the fact that uq ± Ran (P°) in L'^{Q, gW^dx), which follows 
from (11201) . Let f„ be the solutions to 

-V^ ■ cr/j,e(„)V*t;„ = 5',n^(Ao - q)uo, Vnlan = 0. 

From ( |95l) . together with (11221) . (11231) . we see that u„ — — > as 
n oo in L'^{Q, gW^dx). On the other hand, by Prop. 13.11 with Q = 0, 
Vn in L^(f2) as n ^ oo, where Vq is the solution to 

-V'' ■ ctrV^o = ^^/'(Ao - qW Vn\a^ = 0. 

Hence, f o = Mq 7^ is an eigenfunction of Au^rn and belongs in Ran (P°). 
This proves the claim. □ 
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4. Approximate acoustic cloaking 

In this section, we show that, for A ^ spec(^i), it is possible to ap- 
proximate the solutions to the singular, anisotropic acoustic equation 

Aiu-Xu = f, f e L'^{Vl), 

by solutions to certain non-singular, isotropic acoustic equations. Namely, 

Theorem 4.1. Let K G C be a compact set such that K (1 spec{Ai) = 
andXeK. Then, for f G L'^{n), 

hm hm lim(^K,^,, - A)-^ = (A - A)-^, 

it— >1 m— >oo 

in the strong topology of L'^{Q) and weak topology of Hg{Q), uniformly 
with respect to \ ^ K. 

Proof. By Lemma [231 (11 111) , it follows that 

(124) lim(^^,„,,, - A)-V = i-AR^m - A)-V 

in the strong topology of L'^{Q) and weak topology of H^{Q), uniformly 
in A G -ft'. Using Lemma [2.111 we obtain 

(125) lim {Ar,^ - A)-V = {Ar - A)-V 

in the strong topology of H^{Q), uniformly in A G -ftT. Using ([8]) and 
iQ, the convergences in (I124p and (I125P are valid also in Lg{Q) and 

As / G L^{Q) C ^ifl) we have by Lemma [231 

hm(^«-A)-V = Mi-A)-V, 

in the strong topology of Lg{Q) and weak topology of Hg{Q), uniformly 
with respect to A G K. This equation, together with (11241) . (I125p . 
proves the claim. □ 

Let us turn our attention to the convergence of the DN maps. For 
R> 1 and e > 0, the DN map for the acoustic equation, 

(126) -V'-aR,,\/''u + g'/\u-Xgl^\ = 0, u\an = h, 
is the operator 

(127) A^_,_^ .,h^{d, + tiy &)<,,^|aQ. 

Here M^^g^^ the unique solution to fll26p assuming A ^ spec(.4ij,e,m) 
and is the normal (radial) derivative at dfl. 
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Recall that by Lemma [2731 the DN map A^^^ defined in ( H2|) corresponds 
to the boundary measurements for equation flTUl) . In the following we 
assume, for technical simplicity, that b G C^{Q, M^). 

Corollary 4.2. Assume that the magnetic field b is C^-smooth on Vt 
and h G H^^'^(dn), [a, 6] fl spec{Ai) = 0. Then, uniformly with respect 
to \ E [a, b], 

(128) hm hm lim (A^_^,,/i - K\h) = 0, 

where the convergence is strong in the H^^'^(dil) topology. 

Proof. Let G H'^{n), supp(/i^) C {x : 7/3 < |x| < 3}, be an 
extension of h, that is, 

(129) h'^ldQ = h, ||/i''||f/2(Q) < C, if \\h''\\H3/2(^gQ) < 1. 
Then the solution u\^^{h) to fll27p may be represented as 

where 

(130) = -{Ai - \)-^F\h), 

F\h) = -V^- V^h' - Xh' G L^in), supp(F\/i)) c n \ B(]-). 

Here we use the fact that, in the layer fl\B{^), a^^^^ = = 6''^ , g^ = 
g = 1- By Theorem 14. 

lini lim limjv^ (h) - v^(h)\\L2({7/3<\x\<3}) = 0. 

It then follows from [19, Thm. 9.13] that 

lini lim limjv^ (h) - v^{h)\\H^{8/3<\x\<3}) = 0, 

SO that, by trace theorem for Sobolev spaces, we have 

(131) lim lim lim{d^ + iu-b) (u^^^ ^{h) -u^{h)) Ian = 0, 

strongly in H^/'^{dn). □ 

Our next goal is to show the existence of a sequence of non-singular iso- 
tropic mass densities and non-singular bulk moduli, uniformly bounded 
(in n) from above but not from below, with the corresponding acoustic 
operators approximating the singular acoustic operator — A. 
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Theorem 4.3. There exist sequences R{n) 1, m{n) oo, and 
e{n) such that, for any f G L'^{^), A G C, /i G M, A, /i ^ spec{Ai), 

(132) lim(^^(„),^(„),(„)-A)-V = (A-A)-V in LHQ), 

n — >oo 

(134) lim dim{NR^n),m{n),s(n){lJ')) = dim{Ni{fi)). 

n — >oo 

For compact K G C with K fl spec{A) = 0, for n sufficiently large one 
has K n spec{AR[n),m{n)ein)) = o-^d thc limits are uniform in X E K 
and fi e K nR. 

Proof. Let {fp : p = 1, 2, . . . } be a dense set in L^(r2). Let Kn C 
C, = 1, 2, . . . , be a sequence of compact sets so that 

(135) K„nspec(A) = 0, ^nCi^r+i, 

oo 

[jKn = C\ spec(A), dist(spec(^i), K^) > n"^/^. 

n=l 

As L^(r2) C L'^gi^), it fohows from Lemmas 12.81 and \2M that for any 
G Z_|_, there exists a R{n) > 1 such that 

(136) ||(^^(„,) - A)-Vp - (^1 - A)-VpL < K 

dim{Nji(n){fi)) = dim(A^i(/i)), dist(spec (^^j(„)), ir„) > n~^/y2, 

Kin, - Prhl^Ll < ^, K,^,f, - P^fM ^ ^' 

for all p = 1, . . . ,Pn, A G Kn, /i G Kn n M. Here p„ is defined so that 

{/ G L\Q) : <n}c iVi/„({/4^"=J, 

where A^e(5') denotes the e— neighborhoods of S in L^(f2). 

Having chosen R{n), using ([8]) and Lemma [2. Ill and Corollary 12. 12[ we 
choose m{n) such that 

(137) ||(^iJ(„),m(n) - A)"Vp - (^/J(„) - A)"Vpll9 < 

dim(A^H(n),m(n)(/i)) = dim( A^i (/i) ) , 
dist(speC (^R(„),m(n) ),Kn)>n-^/y3, 

for all p = 1, . . . ,pn, A G K^, /X G Kn H M. Having chosen R{n) and 
m{n), one can now use (IE]) and Lemmas 13.31 and [3^ to choose e{n) > 



\'^R{n),m{n),e{n)fp ^R(n),m{n) fpW ^ 



\'^R(n),m{n),e{n)fp ^ifpll^ — ^^ 
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SO that, 

(138) ||(.Ai?(n),m{n),e(n)-A)~Vp-(-4.R(n),m{n)-A)~Vpllg<-; 

dim{NR(^n)Mn),ein)if^)) = dim( A^ij(„)^^(„) (/i) ) , 

dist(spec(^R(„),m{n),e(n)),ifn) > U'^^^/A, 

1 

n 

for all p = 1, . . . , pn, A G -ft'n and /i G fl M; and, 

(139) ||(^i?„,n,e - >^y'\\LHn,gU'd.)^LHn,gU'd.) ^ ^^^^^^ A G 

c/. 0137P and Lemma 13.31 (i). 

Clearly, equations (11361) . (11371) and (I138p imply that 

3 

(140) II (-4.i?(n),m{n),e{n) -A)" /p-(^l-A)" fp\\g<—, 

Th 

dim{Nii(^ri),m(n),e{n)il^)) = dim( A^i (/i) ) , 

3 
n 

for all p = 1, . . . , Pn, A G Kn, yU G -ft'n fl M. 

In particular, (11401) implies (I134p due to UnKn = C \ spec(^i). It 
remains to show that, for / G L^(f2) and a compact set K such that 
K n spec(^i) = 0, one has that {AR[n),m(n),e(n) — ^)~^f converge to 
(^1 — A)~^/ in L'^i^) and that convergence is uniform for all A G -R'. 
Indeed, this will provide equation (I132p . Using the Riesz formula for 
-P_R(n),m(n),e(n)' ^1 ^ ^^^^ ^^^^ ^^^^ pioYe equatiou 01331). 

Let uq G be such that K C Kn for all n > tiq. Furthermore, for 
n > uq, there exists /p(„) with p{n) < p„ such that 

(141) \\f - fpin)\\LHn) 
Clearly, for X e K, 

(142) ||(^,j(„),„.(„),,(„) - A)-V - (^1 - A)-VL 

< ||(-4/?(n),m{n),e(n) - fp(n) " (A " A)"Vp(n)||3 

+ ||(-4i?(n),m(n),e(n) " A)"^(/p(„) - f)\\g + ||(^l - A)"^(/p(„) - /)||g. 

By (I140p . the first term in the right-hand side of (11420 is bounded by 
3/n. To estimate the second term, we use the estimates ([SS]), (I139p . 
and (I14ip . which imply that this term s bounded by 12(2/n)^/^ for 
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n > no. Finally, by (IHD and (|T35ll . 

for A G -ftT. Thus by (11411) . the third term on the right-hand side of 
(11421) is bounded by (S/n)^/^ for n > uq. 

Summarizing, we see that the left-hand side of (I142p tends to 0, as 
n — s> oo, uniformly for X E K, thus proving (11321) . □ 

In the sequel, we will use abbreviated notations 

(143) 0"(n) := '7_R(n),e(n), g(n) '■= gm{n), -^(n) '■= '^R{n),m(n),e(n) 

for the sequences R{n),m{n),e{n) obtained in Theorem 14.31 For sim- 
phcity, denote by A^;^), rather than AR^n),m{n),s{n)y corresponding 
DN map. 

Regarding the convergence of the DN-maps, by the same arguments as 
in proof of Corollary 14. 2[ we obtain 

Corollary 4.4. Assume that the magnetic field b is -smooth in Q. 
Let K C C be compact and such that K fl spec{Ai) = 0. Then, for any 

heH^/^{dn), 

(144) lim Aj„)/i = Al,h, 

where the limit is in H^^'^{dVL), uniformly with respect to \ E K. 

In particular, when b and q\Q\B{i) vanish, A^^^ corresponds to the mea- 
surements on the boundary of a homogeneous ball. Thus (11441) means 
that the isotropic material parameters cr(„) and approximate an 
acoustic invisibility cloak as n — > oo. 



5. Approximate quantum cloaking 



The results of the previous sections can now be used to obtain approxi- 
mate quantum cloaking at a fixed energy, for any potential Q E L°°{Q) 
supported inside the cloaked region -B(l). 

In sequel, we assume that the conductivities aji satisfy fl23l) and that 
also an E C^iU-.R"^). This yields that also aR^,,a(^n) e C^iU-.R"^). 

Let ii^ G M be a given energy level and 

(145) q''{x):=]Q{x) + ^Exi{x), 
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where xi is the indicator function of the ball -B(l). In this section we 
assume that E and Q are such that E is not in the spectrum of the 
operator Ai = Af defined in ( IT9!) using the potential q^. Observe that 
spec(^f ) = spec{Aout) U spec(y4^), where Aout and are defined by 
( HOj) in (14T!) using potential q^. Note that Aout is independent of E and 
Af^ = ^Sin + jE, where Sm is the Schrodinger operator 

S,n = -A + g, V{S,n) = {V2 e H\B{1)) : d,v\9Bil) = 0}. 

Observe that 

(146) E ^ spec(^f ) E ^ spec(5'i„) and E ^ spec{Aout)- 
Then for n large enough, E ^ spec(^(n)). Next, consider the solutions 

Un of 

(147) (-(^(^f V'-a(„)V'' + g^-E)n„ = in fi, u^U = h. 

This equation can be converted to a Schrodinger equation with mag- 
netic potential using the gauge transformation, 

(148) = ali^{x)unix). 
Then 

(149) V'- CT(„) V'(a(;)/Vn) = - VVn + W^nV^n, 

where 

(150) = cr(-;f A (ali^), supp(W^(„)) C 5(2) \ B(l). 

Thus, using the transformation fll48p we see that the acoustic equation 
(11471) for Un is equivalent to the Schrodinger equation for ipn, 

(1/2 \ 
(Tin) 'J 

where for the boundary condition we use ipnlan = Un\dn = h, since 
(T(„) = 1 near dfl. 



Next, define the cloaking potential 

1/2/ \ . 

fjn^W fSE 



1/2/ N 

(152) Kf(x) := l^(„)(x) + ^^(^xi(x)-i5;)+i? 



W^n)ix)+Eil-xiix)){l~^-^l\ 

0{n) \X} 
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where we have used the fact that g^^^^ (x) = 8, o"(„)(a;) = 2 for |x| < 1. 

Thus Vj^ = in -B(l). Clearly, vanishes also near dfl. Then (I15ip 
can be written as 

(153) (-V^- + (x) + Q-E)^r. = 0, Alan = h. 

Now we are ready to prove our main result concerning approximate 
cloaking in quantum mechanics. 

Theorem 5.1. Assume that Q G L°°{Q) is a function supported in 
5(1), b G C^(il,]R^), and E e M are such that E ^ spec{Aout) U 
spec{Sin). Then for any h G H^^^{dQ), 

(154) hm A^.+q/. = Af> ^nH'/\^n). 
Here A^^ , n the DN maps, 

^y„^+Q ■ h ^ d.^nlan; where{-V^- V^+V^+Q-E)^^ = 0, ^n\dn = h, 

while the DN-map A^^ corresponds to the operator Aout with ki = 0, 
see gf . 

Proof. By the hypotheses of the Theorem, it follows that E ^ 
spec((^f), where Ai defined by ( |T9l) with magnetic potential b and 
the potential given by f ll45p . Thus, the Dirichlet problem f ll47p is 
uniquely solvable for large n. As the gauge transformation (11481) 
is the identity map near dfl, we see that d^ipnlan = d^Unlan and 
'0n|eQ = Union- The DN maps for the Schrodinger equation (11531) 
and for equation (11471) thus coincide, and the assertion follows from 
Corollary HH □ 

Note that Theorem 15.11 is of a very different nature than the well- 
known results from the classical theory of spectral convergence, since 
the cloaking potentials do not tend to as — oo. On the contrary 
as seen from the construction of a^^e in Sec. 13.21 and definition (I150p . 
(I152p . supj. iKf (^)l ^ oo as n ^ CX3. Moreover, is of a highly 
oscillatory nature in -8(2) \ B{1) with quasiperiod tending to as n ^ 
oo. 

Theorem 15.11 has two important physical consequences; see ES] for 
further discussion and applications. Consider separately two cases: 

• Suppose that 6 = 0. Since Q is supported in -B(l), the oper- 
ator Aout = —A is the free Schrodinger operator. Then, from 
a physical point of view, the potentials + Q can be con- 
sidered as almost transparent potentials at energy E. Also, the 
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Vj^, which depend on E but are independent of Q, serve as ap- 
proximate invisibility cloaks for two-body scattering in quantum 
mechanics. As all measurement devices have limited precision, 
we can interpret this as saying that, given a specific device, 
one can design, for a given energy level E, a potential to cloak 
an object i.e., an arbitrary potential, from any single-particle 
measurements made at this energy E. 

• Now suppose that Q = 0, while b e C\n,R^), b ^ . We now 
have Aout = — V^^ ■ V^^, and, due to the transformation rules 
for the magnetic potentials, jSi is in general no longer bounded 
near O. Thus, the potentials act as devices which give an 
external observer the illusion that, as —>■ oo, the magnetic 
field is unbounded near O. In particular, for b{x) = Bq x x, 
Bq E M.^, corresponding to a homogeneous magnetic field, the 
illusion Pi{x) has a singularity of the order at O; see [25] 
for details. 



6. DN MAP NEAR EXCEPTIONAL VALUES OF E 

Theorem 15.11 shows that the behavior of A^^irn when E is far from 
spec(^^) and n is large, well approximates the behavior of A^^. This 
situation changes dramatically when E is close to an eigenvalue of the 
cloaked region. 

For simplicity, let us consider the case of an G M which, while 
being in the resolvent set of Aout, is also a simple eigenvalue of Af^ , 
that is, G spec(5'i„), see fll46p . This implies that E^ is a simple 
eigenvalue of Af°. The corresponding eigenfunction then satisfies 
supp C B{1), i.e., vP is a trapped state supported in the cloaked 
region, cf. Lemma [2.31 and (H3|) . 

In the following, let dG(0, 1) be such that dist(E°, spec{Af)\E'^)>d. 

Theorem 6.1. Let Q G L°°{Vt) be a function supported in B{1), 
b G C^(i7,]R'^), and E^ E W be an eigenvalue of Af° corresponding 
to potentials b and defined in Assume that E^ has multi- 

plicity one, and let vP be the corresponding eigenfunction supported in 
5(1). Then 

(i) There is a sequence E(_n), E(n) E^ as n ^ oo, such that E{n) are 
simple Dirichlet eigenvalues of the Schrddinger operators — V**- + 
VinC^+Q- Moreover, the L'^iVt) -normalized eigenf unctions (pt^n) of these 
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Schrodinger operators for the eigenvalues E^n) satisfy, for any p > 2, 

(155) Jim = 0, m C\Tl\B{p)). 

(a) Let h E H^^^(d^l) and ip{^n){h) he the solution to M5^) for some E 
with < 1^ - < d/2. Then 

(156) {d^, + iv ■h)^(^ri){h)\an 

(157) al){h) 

and functions pfn){h) are uniformly hounded, as n oo, in H^f^^dVL)) 
for \\h\\H3/2 <l andEe {E^ - d/A, E^ + d/A). 

Proof. For a given potential Q, the potential defined in (11451) 
depends on E. Thus we start by analyzing how the eigenvalues and 
eigenfunctions of Af and A^^-j change relative to the variation of E. 
Denote hj \{k, E), k = 1,2, . . . , the eigenvalues of Af, numbered in in- 
creasing order and taking multiplicity into account. The dependence of 
these eigenvalues on E then follows from (11461) . Similarly, let X(^n){k, E) 
be the eigenvalues of A^^^ where A^^-^ is the operator of form (I143p with 
q = q^. Observe that 

(158) Af = Af + {q^-q^'), ^f„)=<) + (g^-g^°). 

Using Kato-Rellich formula [33| Thm. VII. 3. 6] and the fact that \\q^ — 
q^Wioc. = ||E - E\, see (11451) . we obtain that 

(159) |A(„)(A;,E)-A(„)(A;,E)| < ^\E - E\. 

In the future we will consider only the value k = such that X{k^,E'^) = 
E^, writing, e.g., \n){E) for A(„)(/c°, £■). Next, consider the spectral 
projectors for p = E^ — d/2 or p = E^ + d/2. Let P^{E) be the Riesz 
projectors for the operators Af and P^^-^{E) be the projectors for A^^y 
They are defined analogously to fl62|) using a contour F C C that sur- 
rounds all of the eigenvalues smaller than p. We can assume that F 
is such a contour that, for n large enough, the distance from F to the 
eigenvalues of the operators Af^ and A^°-^ is more than d/A. Then 

the norm of (-4.^° — z)~^ in LF'{VL, g^^j^dx) is bounded by A/d. Thus, 



= (d^ + iJ' ■b)ip(n)\dn+pfn)ih), 

Jan du 
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assuming that \E — E'^l < d/8, we obtain, using the formula 
in fl62l). the estimate 



(160) \\pi:^^{E) - pjo(^°)iiL^(.j/^..)-^^^(.(^i^^.') ^ ^^'1^ - 

where depends only on the choice of F. 

As is the only eigenvalue of Af° in the interval {E^ — d,E^ + d) 
and has the multiplicity one, it follows from Theorem 14.31 that when 
n is large enough, then A^^-^ has only one eigenvalue E^ in the interval 
{E^ - 3rf/4, E° + 3rf/4) and - < rf/4. Moreover, the eigenvalue 
E^ is simple. 

Let us show that there are -E(n) — ^ such that 

(161) A(„)(i?(„)) = 

Observe that, by Theorem 14.31 A(„)(i?°) — »• as n — > oo. Together 
with f ll59p this implies that for any e > there is n{e) such that for 
n > n{e), 

A(„)(-) : [E' -e,E' + e]^ [E' - e, E' + e]. 

As A(„)(-) is a contraction, see fll59p . we conclude by the Banach fixed 
point theorem that there is a unique E^n) satisfying (11611) . 

Returning to equation fll60p we see that 

(162) Jim ||p,':,(E,„,) - P^.C'") II. = 0. 

Combining this with (11330 (for q = q^°) and embedding fl66|) . we see 

(163) hm P'; (E(„))/ = Pn^°)/ in LliQ) 

Let /i = E+d/2 and z/ = E—djl. For large n, the operators P^'^J^^ (£'(„) ) — 

PI^^^{E(^n)) are the orthoprojectors, in LF'{Vt^ g\j^dx) onto the eigenspace 

of corresponding to the eigenvalue -E(n). Using (11620 we see that 
for n large enough, the eigenvalue E(^n) has multiplicity one. Denote 
by U{n) the eigenfunction corresponding to normalized in L'^gi^t). 
Now M° is supported in P(l) and thus 
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Using (11631) we see that 

in where a„ = |a„|e*"" and |a„| — > 1 as n — oo. 

Denoting = e*""M(„), we see that 

(164) hm uin) = M° in Ll{n). 
Since M°(a;) = 0, |x| > 1, this imphes 

(165) Jirn \\u(^n)\n\B{2)\\L^n\B{2)) = 0. 

Observe that (T(„) = 70, = 1 and V^^^ = in Q \ B{2). Thus it 
follows from (11651) that the functions satisfy 

-V^- V''m(„) = in \ B{2), U(n)\dn = 0, 

where the right side E(^n)U{n) ^ in L'^{Q\ B{2)) as n 00. Since 
6 G C^{Q \ -8(2)), standard elliptic regularity results [19] imply that 

(166) lim ||M(„)||c>i(n\B(p)) = 0, 

n — 5- 00 

for any 2 < p < 3. Using the transformation (11481) to define (pnix) = 
(7^^^^ {x)un{x) , we see that 

(-V''- V' + Vn'"' + Q)0„ = (pnlan = 0. 

This proves that 0„ is an eigenfunction of the Schrodinger operator 
—V^- + Ki + Q for the eigenvalue -E(n)- Moreover, since (T(„) = 2 
in 5(1), it follows from (ED and ffTO that 

liminf ||0„||i2(m > a/2 liminf 1^^2(5(1)) > V2\\u^\\ 1^(^3(1)) = J- 

n — ^00 n — *-oo Z 

This inequahty and equation (11661) together imply (I155P . Thus we have 
proven (i). 

Next, consider (ii). We start with the boundary- value problem for the 
acoustic equation 

V"- a^V^ + Vf + q^ + E) v^ih) = 0, v^ih)\sn = h. 

Clearly 

(167) v^{h) = a^{h)u(^n)+w^{h), with (wf(/i), u^n)) ^2^^^gi/^'2^dx) = 0' 
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where, using the notations introduced in (11291) . (11301) . we have 
(168) af (/i) = (/^^ + ^^^^ {F'^W, , 

wl{K) = (^(„,) - E)-' (^F^(/i) - (F^(/i), W(„))^,(^v.^^.) . 

Since F-^{h) — {F^{h), U(^n)) M(n) and U(^n) are orthogonal in L^(fi, g^^^^dx), 
and dist(ii^, spec(^(„)) \ {-E(„)}) > (i/2, it follows from (11291) that, for 
n large enough, 

where C is independent of G {E^ — d/4, E^ + d/4) and h satisfying 
||/?.||j:^3/2(gj^) < 1. Note that in fl\B{2) the function w^{h) satisfies the 
equation 

w^{h)\9n = 0. 

Thus, by boundary elliptic regularity, |l9j, Thm. 9.13], 

\\w^{h)\\m{n\B(p)) < Cp, 2<p<3. 

This inequality, (I129P and (11301) together imply that pfn)W = {(^u + 
iu ■ h)w^{h) satisfies 



(169) \\P{n)ih)\\m/2(^s^) < Co, for \\h\\HV\dn) < 1, 

if n is large and \E — -E(ri,)| < d/2. Finally, integration by parts shows 
that 



gf (/.) = ^ \ I h?^dS= , \ [ h?p^dS. 

The desired equation (I156p follows from the above equation together 
with (I167p . (I169P and (I157p . if we take into the account the relation 
(fT48D between ^f^){h) and v^{h). □ 

Remark. Theorem 16.11 means that, away from the cloaking structure, 
the eigenfunctions converge to zero as n — cx), i.e., the rep- 
resent almost trapped modes, effectively vanishing near dVt. Physically 
speaking, we can say that if E^ is an eigenvalue of the Schrodinger op- 
erator with Neumann boundary condition in the cloaked region -B(l), 
and we connect the interior and the exterior via the cloaking potential 
in the layer B{R(^n)) \ -B(l), with — ^> 1 as n — > oo, a particle 
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under the influence of tlie combined potential on -8(3) is still largely 
confined (modulo standard tunneling) to the interior region, with a 
slight shift of the energy of the eigenmode from to 

Moreover, for energies E close to the values the presence of the 
cloaked region is very clearly seen in the boundary measurements of 
the DN map, so that the invisibility effect is compromised. On the 
other hand, at energies which are away from the -E(n), the DN map for 
the potential Vj-f^ + Q well approximates A^^, and thus the potential 
Q is approximately cloaked. 



7. Numerical results 




Figure 1. Solutions of acoustic equations. Left: 
Solution of the scattering problem (11701) restricted 
to ball -8(3), when a plane wave scatters from an ap- 
proximate cloak in the case fll75p . i.e., when fc^ is away 
from the exceptional values -B(n) . Right: Almost trapped 
eigenfunction u of the acoustic operator (11721) . with 
Dirichlet boundary condition, /i = 0, in the case (11761) . 
i.e., when k"^ is equal to the exceptional value E(^n)- 

Next we consider scattering problems for the Helmholtz and Schrodinger 
equations in the case when the magnetic potential vanishes, i.e., b = 0. 

The scattering problem for Helmholtz equation is 

(170) (V-a(„)V + P(l + a(x))(7f/)V* = in M^ 

u'"\x,k) =u'''{x,k) + u'%x,k) 
and for the Schrodinger equation 

(171) (-V- V + v;^ + g - E)^*°* = in m.\ 

^*°*(x, E) = V^"(x, E) + iP'^^ix, E), 
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where E > 0, k = E^/^, the incident fields are u"'{x, k) = il)"\x, E) — 
^tkui-x^ _ 2 and the scattered fields satisfy the radiation condition 

d d 
lim r(- ik)u'"'(x, k) = lim r(- iE^^'^)ib^''(x, E) = 0, r = \x\. 

r-*oo or r-^oo Or 

In following we consider a{x) that corresponds to a real bounded po- 
tential Q supported in B{1), that is, a{x) = —{E~^Q[x) + 3)/4. We 
assume that Q{x) is such that 1 + a{x) > cq > 0. 

We consider also the solutions of the boundary value problems in Q 
(note that we denote these solutions by u and ip, without using super- 
scripts) , 

(172) (v-a^n)V + k^{l + a{x))gll^^u = Q in Q, u\9n = h 
and for the Schrodinger equation 

(173) {-V-V + V^ + Q-E)il;^0 in Q, u\9a ^ h. 

The solutions of these scattering and boundary value problems are 
related through a gauge transformation, 

(174) i;'"\x) = ali^{xy"\x), i^ix) = ali^){x)u{x). 

The computations are made without reference to physical units; for 
simplicity, we use E — 2. The cloak corresponds to the parameter 
R = 1.005 and inside the cloak we have located a spherically symmetric 
potential; 

Q{x) = QinX[o,R]{\A)^ that is, a{x) = -(5i„(4£;)-^X[o,ii](k|) - 3/4. 
To illustrate the approximate cloaking, we used 

(175) Qin = 1, 
and to obtain an almost trapped state, 

(176) = -2.576. 

In our numerical solution we have approximated a(^n) by a picccwisc 
constant function consisting of 30 layers in the region R < r < 2. The 
values of the conductivity in these layers are chosen as in the above 
sections. This corresponds to the case when the cloaking potential 
is a weighted sum of delta functions, and their derivatives, on spheres. 

In the numerical solution of the problem, we represent the solution 
and u in terms of spherical harmonics and Bessel functions up 
to order N = 7 in each layer where the cloaking conductivity is con- 
stant. The transmission condition on the boundaries of these layers are 
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Figure 2. Radial profiles of solutions of acoustic 
and Schrodinger equations. Left: The solution 
of the scattering problem (11701) on the line segment L = 
{(a;,0, 0) : x G [0,3]} in the case of Fig. l(left), where 
is far from the exceptional values -E{n), is shown with 
the blue curve. Also, the eigenfunction u of acoustic 
operator fll73p on the line segment L in the case of Fig. 
1 (right), where k"^ is the exceptional value -E(n), is shown 
with the red curve. Right: The solutions and ip of 
the Schrodinger equations (11711) and (I173P on the line 
segment L, obtained from the solutions on the left via 
gauge transformations (11741) . 



solved numerically by solving linear equations. After this we compute 
the solution -0*°* and of the Schrodinger equation using the gauge 
transformation (11741) . 
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